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ABSTRACT 


The  decision  problem  of  the  allocation  of  available  weight  and 
space  in  an  air-dropped  anti-submarine  torpedo  to  fuel^  explosive^,  and 
to  attack  speed  capability  above  some  preset  minimum  is  studied  for  the 
case  in  which  alternatives  under  consideration  do  not  differ  with  respect 
to  their  effect  on  weapon  reliability  and  to  their  effect  on  enemy 
countermeasures  capability  *  Analytical  models  are  developed  which 
relate  (l)  probability  of  acquisition  to  endurance  for  the  case  of 
circling  search  at  constant  depths  (2)  probability  of  hit  given  acqui¬ 
sition  to  endurance  and  attack  speedy  (3)  probability  of  hit  on  a 
close-in  attack  run  to  attack  speedy  and  (k)  probability  of  kill  damage 
given  a  hit  to  amount  of  explosive  usedo 
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CHAPTER  I 


INTRODUCTION 

The  Problem 

An  air-dropped  anti-submarine  torpedo  may  be  characterized  by  a 
long  list  of  parameters o  Specification  of  parameter  values  in  a 
specific  weapon  is  made  such  that  the  weapon  may  achieve  a  maximum 
of  effectiveness  in  an  actual  tactical  situation  in  wartime 0  The 
state  of  the  art  in  weapon  design  will  determine  the  values  of  many 
of  the  weapon  parameters*  For  example*  the  transducer  parameters 
may  be  prespecified  by  the  status  of  some  continuing  torpedo  sonar 
development  program.  Cost  considerations  may  also  be  of  great 
importance  in  the  evaluation  of  alternative  proposals  for  the 
design  parameter  values.  For  example*  a  choice  between  two  propul- 
sion  systems  may  hinge  on  a  difference  in  development  costs  necessary 
to  produce  a  reliable  weapon.  However*  for  torpedoes  and  more 
especially  for  air-dropped  torpedoes*  the  restrictions  of  weight 
and  of  physical  dimensions  (or  volume)  resulting  from  the  varied 
mix  of  capabilities  of  the  prospective  weapon  delivery  systems  are 
also  very  significant  in  the  weapon  design. 

Consider  the  restraints  of  weight  and  volume.  One  of  the  design 
parameters  which  has  a  great  effect  on  weapon  weight  and  volume  is 
the  maximum  operating  depth.  An  increase  of  a  few  hundred  feet  in 
the  operating  depth  requirement  will  require  a  more  substantial 
structure.  This  in  turn  may  result  in  a  decreased  weight  and  volume 
to  be  made  available  for  other  weapon  components.  However*  intelli¬ 
gence  reports  or  estimates  of  the  depth  capabilities  of  prospective 


submarine  targets  in  the  time  period  for  which  the  weapon  is  to  be 
designed  may  prompt  the  decision  maker  to  fix  the  design  operating 
depth  as  an  input  to  the  design  problem e  Similarly  the  advance  speei 
fication  of  a  minimum  attack  speed  capability  may  be  prompted  by 
intelligence  estimates  and  will  certainly  act  to  place  an  upper  limit 
on  weight  and  volume  available  for  other  weapon  components*,  Suppose 
that  a  study  of  specifications  from  higher  authority^  engineering 
capabilities  to  produce  the  various  weapon  components^  cost  considera 
tions  or  some  appropriate  optimization  procedures  have  resulted  in 
specification  of  values  for  all  weapon  design  parameters  except 
1*  Amount  of  fuel 

2.  Amount  of  explosive 

3.  Maximum  speed  (above  some  preset  minimum  value  )c 

Further,  suppose  the  restraints  on  weight  and  volume  have  been  deter- 
mined  such  that  if  the  amount  of  fuel  and  explosive  and  the  maximum 
speed  capability  above  some  specified  minimum  are  considered  as  a 
weapon  payload,  then  the  available  payload  weight  and  volume  are 
provided  to  the  decision  maker 0  The  decision  problem  of  allocation 
of  available  payload  weight  and  volume  to  amount  of  fuelp  amount  of 
explosive,  and  to  maximum  speed  capability  above  some  specified 
minimum  is  the  problem  considered  in  this  paper 0 
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The  Approach 

An  initial  approach  to  the  allocation  problem  is  to  set  certain 
minimum  limits  on  the  parameter  values .  The  minimum  value  of  maximum 
(attack)  speed  has  been  specified 0  There  must  be  sufficient  fuel  to 
provide  some  minimum  probability  of  target  acquisition0  (in  the  case 
of  continuous  circular  search  at  a  constant  preset  depth  as  is  con¬ 
sidered  in  this  paper  a  minimum  probability  of  acquisition  might 
result  from  a  dive  to  search  depth  and  completion  of  one  complete 
search  turn.)  Then  there  must  be  sufficient  fuel  for  at  least  one 
attack  under  some  assumed  conditions  of  acquisition  range  and  target 
course  and  speed.  Finally  there  must  be  at  least  enough  explosive 
to  result  in  kill  damage  if  a  hit  occurs  at  the  most  favorable  depth 
and  point  of  contact  on  the  target  surface.  Having  determined  these 
allocation  minimums  there  remains  the  problem  of  allocating  the 
balance  of  available  payload  weight  and  volume. 

A  second  step  may  determine  a  maximum  value  of  attack  speed. 

This  maximum  may  be  determined  from  a  study  including  the  probable 
operating  depths  the  characteristics  of  the  guidance  system^,  and  the 
generated  self-noise  level  and  its  effect  on  the  ability  of  the  sonar 
to  maintain  contact  with  the  target  on  the  attack  run. 

What  is  the  next  step?  How  does  the  decision  maker  decide  between 
enough  additional  fuel  for  one  reattack  and  enough  additional  explosive 
to  produce  kill  damage  over  a  given  additional  element  of  target 
surface?  What  measure  of  effectiveness  is  used? 

This  paper  proposes  a  procedure  which  may  provide  the  decision 
maker  with  a  quantitative  analysis  or  comparison  of  the  tradeoffs 
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involved  in  alternatives  of  allocation  among  fuel,  explosive,  and 
additional  attack  speed  capability 0  The  tradeoffs  are  measured  in 
terms  of  probability  of  kill  as  a  measure  of  effectiveness *  Probabil¬ 
ity  of  kill  is  dependent  upon  weapon  reliability  and  upon  the 
capability  of  the  enemy  to  employ  countermeasures  (in  addition  to 
maneuvering)*  It  is  assumed  in  this  paper  that  weapon  reliability 
and  enemy  capability  to  employ  countermeasures  (in  addition  to  maneu= 
vering)  are  both  independent  of  the  amount  of  fuel  and  explosive  and 
of  the  weapon  attack  speed  within  the  range  of  values  between  the 
preset  minimum  and  some  determinable  maximum*  This  enables  removal 
of  reliability  and  enemy  countermeasures  capability  (less  maneuvering) 
from  the  problem  and  probability  of  kill  may  then  be  factored  into 
three  independent  probabilities j  those  of  probability  of  acquisition, 
probability  of  hit  given  acquisition,  and  probability  of  kill  given 
hit*  An  optimization  is  made  of  probability  of  kill  as  a  function 
of  the  fuel,  explosive,  and  additional  attack  speed  capability  with 
weight  and  volume  restraints*  This  does  not  provide  a  number  for 
use  as  an  operational  value  of  kill  probability*  However  the 
dynamic  programming  methods  used  provide  a  solution  which  specifies 
the  optimal  allocation  of  weight  and  volume* 

It  is  emphasized  that  a  quantitative  comparison  of  the  tradeoffs 
involved  in  the  allocation  process  may  be  obtained  only  by  some 
common  analysis  of  the  three  study  variables  as  they  act  upon  the 
events  of  acquisition,  hit,  and  kill* 
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The  dynamic  programming  optimization  model  developed  generates 
an  input  requirement  for  (l)  probability  of  acquisition  as  a  function 
of  available  fuel^  (2)  probability  of  hit  as  a  function  of  available 
fuel  and  attack  speed  and  (3)  probability  of  kill  damage  as  a  function 
of  explosive  weight.  These  requirements  may  be  satisfied  in  a  number 
of  ways  including 

1.  analytical  derivations  in  terms  of  obtainable  or  estimable 
inputs 

2.  simulation  (war  gaming) 

30  extrapolation  of  empirical  data  from  operational  or 
development  tests 
4.  a  combination  of  the  above 

The  model  user  may  choose  any  of  these  methods  so  long  as  the 
resulting  inputs  are  compatible  with  the  allocation  model .  However^ 
any  of  the  methods  used  will  require  an  assumption  about  the  tactical 
situation  with  respect  to  which  the  optimization  is  to  be  made c 
Consider  the  actual  tactical  situation  in  which  the  weapon  will  be 
used.  Assume  that  methods  exist  to  determine  the  optimum  values 
of  all  weapon  parameters^  given  that  the  weapon  will  be  used  in  a 
completely  specified  tactical  situation.  The  optimum  parameter 
values  determined  will  be  different  for  each  different  tactical 
situation  considered.  The  space  of  different  tactical  situations 
is  infinite.  However  a  single  set  of  parameter  values  must  be  spe° 
cified  for  the  weapon^  ultimately  to  be  produced  in  quantity.  By 
what  criteria  is  the  specific  tactical  situation  to  be  selected  for 
use  in  the  optimization  of  the  weapon  parameters?  The  probability 


5 


that  the  weapon  will  actually  be  used  in  the  tactical  situation  selected 
is  zero.  The  following  two  criteria  will  be  used  in  this  paper  in  the 
analytical  derivation  of  the  three  required  probabilities M 
Criterion  1;  Select  an  unfavorable  tactical  situation  which  is  not 
extreme.  This  may  be  thought  of  as  a  modified  fmaximinl  criterion . 

A  maximin  criterion  implies  optimization  of  the  weapon  parameters  to 
maximize  effectiveness  for  that  tactical  situation  which  offers 
minimum  probability  of  weapon  success.  The  extreme  situation  is 
modified  in  this  paper  because  the  low  probability  of  kill  resulting 
from  assumption  of  an  unfavorable  situation  in  the  extreme  would 
undoubtedly  distort  the  relationships  between  the  study  variables  which 
would  apply  in  a  greater  percentage  of  cases.  In  general  a  weapon 
which  is  optimized  for  an  unfavorable  situation  will  be  able  to  cope 
with  a  more  favorable  situation  for  which  its  parameter  values  were 
not  optimized. 

Criterion  2:  Select  a  tactical  situation  which  is  sufficiently 
simplified  as  to  be  within  the  scope  of  this  paper.  Progress  in 
work  of  this  kind  is  more  often  made  by  starting  with  something 
relatively  simple  and  adding  sophistication  after  experience  is 
gained  with  the  model. 
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Relation  of  this  work  to  previous  work 

The  allocation  problem  has  been  studied  before  by  R,  Ce  Brumfield 
in  references  1  and  2  in  which  algebraic  expressions  were  derived 
relating  total  weight  and  volume  of  a  weapon  of  arbitrary  speed  and 
range  to  performance  characteristics  anti  component  characteristics. 
This  enabled  determination  of  optimum  running  speed  and  range  for 
stern  chase  vehicles  under  conditions  of  minimum  propulsion  system 
weight,  minimum  total  weight  of  propulsion  system  and  associated 
shell,  and  minimum  propulsion  system  volume 0  The  approach  in  this 
paper  is  essentially  different  from  that  of  references  1  and  2 
although  the  problem  is  much  the  same*  The  previous  work  optimized 
total  weight  and  volume  and  attack  speed  for  the  event  of  a  single 
attack  run,  given  explosive  weight,  a  length  of  run  in  the  acquisition 
phase,  and  acquisition  range.  This  paper  optimizes  an  allocation  of 
payload  weight  and  volume  to  fuel,  explosive,  and  attack  speed  for 
the  events  of  acquisition,  hit  (attack  and  reattacks)  and  kill 
given  the  maximum  allowable  payload. 

A  great  deal  of  work  has  been  done  in  the  study  of  various 
phases  of  the  torpedo  problem.  References  3  and  k  form  the  background 
on  which  the  development  in  Chapter  III  is  based.  Reference  k 
develops  a  continuous  volume  scan  rate  for  a  single  circular  search 
turn.  In  this  paper,  in  the  first  three  sections  of  Chapter  III, 
a  volume  scan  for  a  single  ping  on  the  first  search  turn  is  developed. 
This  is  a  modification  of  the  continuous  scan  rate  to  the  discrete 
case  where  the  discrete  interval  is  a  single  ping  cycle.  However,  the 
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development  in  this  paper  is  continued  to  obtain  a  probability  of 
acquisition  not  only  for  a  single  circular  search  turn  but  for  con° 
tinuous  turns  until  fuel  exhaustion  *  This  further  development  is 
made  possible  because  a  tactical  situation  is  postulated  in  accordance 
with  criteria  one  and  twoe  The  postulation  is  Chapter  III  limits 
the  usefulness  of  the  specifics  of  the  development  to  the  search 
pattern  considered  but  is  required  in  order  to  obtain  the  probability 
of  acquisition  input  to  the  allocation  model 0  The  logic  of  the 
development  may  be  modified  as  appropriate  to  provide  for  other 
postulated  tactical  situations  in  accordance  with  requirements  to 
study  a  particular  search  pattern 0  The  study  of  search  turns  subse~ 
quent  to  the  first  and  until  fuel  exhaustion  is  necessary  to  provide 
probability  of  acquisition  as  a  function  of  fuel  weight  throughout 
the  range  of  the  argument * 

The  development  in  Chapter  IV  of  probability  of  a  hit  in  any  given 
pass  at  the  target  has  no  known  precedent  in  the  literature  and  is 
presented  as  an  alternative  way  of  looking  at  the  hit  probability 
problem  which  may  be  of  interest « 

It  is  emphasized  that  Chapters  III,  IV,  and  V  are  samples  of 
possible  methods  of  attacking  the  problem  of  finding  inputs  to  the 
allocation  model  of  Chapter  I  which  are  specifically  designed  to 
permit  examination  of  the  tradeoffs  between  alternative  allocations 
of  payload  weight  and  volume  among  the  three  study  variables 0 
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Structure  of  the  paper 

Chapter  I  contains  the  development  of  a  mathematical  model  of 
the  torpedo  problem  in  a  form  such  that  the  dynamic  programming 
techniques  of  Chapter  II  of  reference  5  may  be  applied  to  solve  the 
allocation  decision  problem.  Chapters  III^  IV^  and  V  contain  the 
development  of  the  three  required  inputs  to  the  allocation  model  of 
Chapter  11%  probability  of  acquisition^,  probability  of  hit  given 
acquisition  and  probability  of  kill  given  hit.  Appendix  A  is  a 
description  of  the  computer  flow  chart  notation  to  be  used  in  the 
flow  charts  of  the  other  appendices.  Appendix  B  is  the  computer 
flow  chart  of  the  dynamic  programming  allocation  model  of  Chapter  I0 
Appendices  C  1  through  C  8  contain  the  mathematical  details  of  the 
developments  in  sections  I  through  VIII  of  Chapter  III  plus  the 
corresponding  computer  flow  charts.  This  mathematical  supplement 
to  Chapter  in  is  included  in  appendices  to  enable  emphasis  of  the 
logic  flow  in  Chapter  III  with  a  minimum  of  mathematical  detail. 
Appendix  D  contains  the  computer  flow  chart  for  Chapter  IV. 
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Summary 

A  mathematical  model  is  provided  which  enables  optimization  of 
the  decision  problem  of  the  allocation  of  available  weight  and  space 
in  an  air-dropped  anti-submarine  torpedo  among  three  study  variables® 
amount  of  fuel*  amount  of  explosive  and  attack  speed  capability  above 
some  preset  minimum.  The  model  is  restricted  to  the  study  of  alterna¬ 
tive  allocations  which  do  not  differ  with  respect  to  their  effect  on 
weapon  reliability  and  to  their  effect  on  enemy  countermeasures 
capability  (other  than  maneuvering  capability.) 

Mathematical  models  are  developed  which  relate  (l)  probability 
of  acquisition  to  endurance  for  the  case  of  circling  search  at 
constant  depth*  (2)  probability  of  hit  given  acquisition  to  endurance 
and  attack  speed*  (3)  probability  of  hit  on  a  close-in  attack  run  to 
attack  speed*  and  (U)  probability  of  kill  damage  given  a  hit  to  amount 
of  explosive  used.  These  four  models  are  developed  as  specific 
examples  of  inputs  needed  to  be  compatible  with  the  requirements  of  the 
allocation  model.  Computer  flow  charts  are  provided  for  each  model 
except  the  last  to  facilitate  programming  in  any  desired  computer 
language  and  also  as  an  alternative  method  of  study  of  the  development 
of  the  model  logic. 
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CHAPTER  II 


THE  ALLOCATION  MODEL 

Consider  the  design  of  an  air-dropped  anti-submarine  torpedo  such 
that  the  amounts  of  fuel  and  explosive  to  be  included  and  the  capabi¬ 
lity  of  attack  speed  in  excess  of  some  predetermined  minimum  value 
may  be  considered  together  as  a  payload  for  which  there  is  available 
a  predetermined  maximum  of  weight  and  volume 0  Consider  the  employ¬ 
ment  of  the  torpedo  from  the  time  it  splashes  into  the  water  until 
the  time  it  explodes  in  the  proximity  of  a  submarine  target  as  an 
ordered  sequence  of  the  following  three  events  each  of  which  must 
occur  in  succession  for  accomplishment  of  the  torpedo  missions 

1.  Acquisition  after  a  search  run  following  some  preset  pattern 

2.  Hit  after  an  attack  (and  if  necessary  a  series  of  reattacks) 

3.  Kill  damage  resulting  from  explosion  ot  the  weapon  in  the 
proximity  ot  the  target 

Divide  the  study  variable,  amount  ot  fuel,  into  two  separate  and 
independent  study  variables: 

1.  Amount  ot  fuel  to  be  used  in  tne  event  of  acquisition  (event  l) 

2.  Amount  of  fuel  to  be  used  in  the  event  of  hit  (event  2) 

Restrict  the  study  variable,  attack  speed,  to  a  range  of  values 
such  that 

1.  The  lower  limit  of  the  range  is  some  preset  minimum  value  which 
is  greater  than  the  given  value  of  weapon  search  speed 0 

2.  As  the  attack  speed  is  allowed  to  vary  over  its  range,  the  fuel 
expenditure  rate  at  attack  speed  and  at  any  depth  under  considera- 
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tion  may  also  vary.  However  as  the  propulsion  system  is  modified 
to  allow  for  changes  in  the  value  of  its  maximum  (attack)  speed 
capability,  the  fuel  expenditure  rate  at  search  speed  and  depth 
remains  approximately  constant. 

Make  the  following  assumptions? 

Assumption  1:  Weapon  reliability  and  target  capability  to  employ 
countermeasures  other  than  maneuvering  are  both  independent  of  changes 
(within  the  ranges  to  be  considered)  in  the  values  assigned  to  the 
four  study  variables: 

1.  Amount  of  fuel  for  event  1 

2.  Amount  of  fuel  for  event  2 

3.  Amount  of  explosive 

4.  Weapon  attack  speed 

Assumption  2:  The  weapon  is  100$  reliable. 

Assumption  3:  The  enemy  target  submarine  has  no  countermeasures 
capability  except  that  of  maneuvering. 

Note:  Assumptions  3  and  4  are  made  in  order  to  justify  the  relation¬ 
ship  between  the  study  variables  and  the  measure  of  effectiveness  to 
be  introduced  later.  The  fact  that  they  are  intuitively  incorrect 
does  not  detract  from  the  validity  of  the  optimization  in  the  model 
because  of  assumption  1. 

Assumption  4:  The  three  events  of  acquisition,  hit  and  kill  are 
independent  events  in  a  probabilistic  sense.  That  is,  if 
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P  [xj  =  Probability  that  event  x  occurs 
=  P  (target  acquisition  (event  1 )] 

Pjj  =  P  [target  is  hit  (event  2)  given  that  event  1  has  occurred  j 
Pp  =  P  (kill  damage  is  inflicted  on  the  target  by  explosion  of  the 
weapon  warhead  (event  3)  given  that  event  2  has  occurred 


If  PR  =  P  [target  kill] 

then  P  =  P  • P  • P 
nen  rk  “  rA  rH  E 

Define : 

W  =  available  payload  weight  (lb) 

V  =  available  payload  volume  (ft^) 

=  fuel  weight  for  acquisition,  event  1  (lb) 
Wg  =  fuel  weight  for  hit,  event  2  (lb) 


w^'s  weight  for  modifications  to  the  propulsion  system  to  provide 
a  maximum  (attack)  speed  capability  in  excess  of  some  preset 
minimum  value  (lb) 

=  preset  minimum  value  of  attack  speed  (yd/sec) 
w^  =  explosive  weight  (lb) 

o 

v^  =  fuel  volume  for  acquisition,  event  1  (ft3) 
v^  =  fuel  volume  for  hit,  event  2  (ft^) 

v^1 =  volume  for  modifications  to  the  propulsion  system  to  provide 
a  maximum  (attack)  speed  capability  in  excess  of  Sq  (ftJ) 

O 

v^  =  explosive  volume  (ftJ) 

w^  =  payload  weight  for  hit,  event  2  (lb) 

i  it 

w2  =  w2  +  w2 
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v2  =  Payload  volume  for  hit,  event  2  (ft*) 

V2  ‘  V2  +  v2 

A  w  =  increment  of  weight  of  the  any  of  the  weight  variables^  all  of 
which  are  varied  in  increments  of  equal  size  (lb) 
gi(wi)  =  the  relationship  between  the  i  th  weight  variable  and 

I  V  0 

the  i  th  volume  variable^  i  «  1,  2,  2  ,  3 

t  t  t 

g*^)  =  v£;  i  »  1,  2,  2  ,  3j  gL(w)  =  g^(w) 

f 1 C wi )  =  natural  logarithm  (in)  of  P^  expressed  as  a  function  of 

WV  )  =  In  P„  expressed  as  a  function  of  w*  and  w“  1 
c.  c.  d  H  d.  2 

f,(w  )  =  In  P  expressed  as  a  function  of  w, 
j  J  r  3 

3 

F  =  In  P  =  XI  fi 
i  =  1  1 


I  _i 

w  ,  w  =  lower  and  upper  limits  respectively  on  the  allowable  range 
of  fuel  weight  (lb) 

i  i  _ i  t 

W  >  W  =  lower  and  upper  limits  respectively  on  the  allowable  range 
of  weight  for  modifications  to  the  propulsion  system  to 
provide  a  maximum  speed  capability  in  excess  of  (lb) 


*  i  i  til 

W  W 

—  9  =  lower  and  upper  limits  respectively  on  the  allowable 

range  of  explosive  weight  (lb) 

The  concept  ot  the  lower  and  upper  limits  on  the  ranges  of  the 
study  variables  is  as  follows*  It  will  be  obvious  to  the  design 
engineer  that  very  small  values  or  either  fuel  or  the  explosive 
weights  will  result  in  a  negligible  probability  of  success*  Thus 
there  is  no  need  to  examine  the  vicinity  of  the  extreme  values  of 
zero  and  W .  However  the  values  of  the  limits  may  be  chosen  as 
close  to  zero  and  W  as  is  desired*  In  the  case  of  W° 0 9  the  value 
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chosen  may  well  be  zero  as  this  implies  that  the  considered  minimum  of 

8 

attack  speed  is  to  be  S q,  the  preset  minimum.  In  the  case  of  W  the 
value  chosen  will  probably  be  well  under  W  as  the  maximum  of  the  range 
of  values  of  attack  speed  is  limited  by  considerations  of  resultant 
guidance  problems,  excessive  self-noise  levels,  and  so  forth .  The 
hit  event,  event  2,  contains  two  study  variables  which  must  be  con¬ 
sidered  together.  For  events  1  and  3  only  one  study  variable  is  con¬ 
sidered  per  event  (or  programming  stage).  Thus  for  each  value  of 
fuel  weight  (w^)  there  is  associated  a  single  value  of  acquisition 
probability  P  (w1  )  and  a  single  value  of  volume  v,  *  g^Wp). 

Similarly  for  values  of  explosive  weight  w^« 

However 

W2  =  W2  +  W2' 

T2  =  T2  *  t2*  =  82(w2)  *  82*  (w2  5 

This  requires  that  for  each  value  of  weight  (w^)  there  be  an 

associated  set  of  values  of  hit  probability  ^(W2^  an<^  volume 

v  .  one  value  for  each  combination  of  values  of  w'  and  w'*  which 
2 7  d  d 

♦  tt  j  ®  8 

sum  to  the  considered  value  of  w^.  A  suboptimization  of  w^  and  w^ 
to  provide  a  single  optimal  value  of  )  for  each  w^  as  an  input 

to  the  allocation  model  would  neglect  the  associated  volume 
restraints.  Therefore  index  the  weight  w^  by  k(k  -  1, K)  over 
its  range 

min  £w'  ,  W"”^  ~v2  “mijl  ^W‘  +  ’  »  w  | 
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and  define 


w  k  =  k  the  value  of  weight  to  be  distributed  between  w^  and 


w“;  k  =  I,-.-,  K  (lb) 

W2jk  =  J  f*1  value  W2  assocfatec^  with  j  -  li,‘,°,'<>  (lb) 

W2jk  =  J  value  of  w^'  associated  with  w^j  j  s  (1^) 

P  k  =  J  th  value  of  associated  with  W^j  j  s  1#’°% 
v^a  =  j  th  value  of  v^  associated  with  w^j  j  -  1*  "*>  JkS 


2  jk  J  2 

k  =  1,  •••,  K 


V2jk  =  j  t*ie  va^ue  °f  v2  associated  with  w^j  j  =  1 


v  =  v*  .  +  v'\,  =  j  the  value  of  v.  associated  with 
2 jk  2 jk  2 jk  J  2 

w2ki  J  =  h  Jk  (ft3) 


The  available  payload  W  and  the  upper  and  lower  weight  limits  or 
the  weight  increment  A  w  are  so  chosen  that  all  the  weights  are 
divisible  by  Aw.  Then  there  will  be 


W*  W» 
Aw  ”  Aw 


considered  values  of  w^  denoted  w^j  r  =  1,  R^, 

there  will  be 


min  {  W»  +  W"  ,  W J-  min  £  W» ,  W«  *j- 


Aw 


Aw 


considered  values  of  w^  denoted  w^J  k  =*  1, 


-  K 

K 
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and  there  will  be 


W"  W« " 


Aw 


Aw 


=  S 


considered  values  of  denoted  w_  s  s  =  ls  *  * ° 9  S 

3  3s"  * 

The  functions  g^w^,  g^Cw^),  g^'  (w^'  ),  g^w^)*  f^Wj),  and 
f^(w^)  are  assumed  to  be  non»decreasing 0  That  is^  an  increase  in 
the  argument  will  not  decrease  the  value  of  the  function 0  In  the 
case  of  fg(wgj  w^1 )  it  is  assumed  that  any  increase  in  either  of 
the  arguments  while  maintaining  the  other  fixed  will  not  decrease 
the  value  of  the  function. 

The  problem  is  to  maximize  F  =  In 
subject  to 

<  W 

Vi  *  v2  +  v3  <V 

Consider  the  new  problem  to  maximize 

p»  =  F  -  X  j^g1(w1)  +  g£(w£)  +  g£’  (w£>)  +  g3(w3)] 


subject  to 


+  w2  +  <W 


gi^wi^  1  =  2»  2  9  3 

where  X  is  a  non-negative  real  number  (Lagrange-multiplier  ) , 


Define  the  recurrence  relations! 


h1(w)  =  fx(W)  -  X  g2(W)'  if  W  =  ^  T  =  1,  •••, 


R 


=  -  co 


otherwise 
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h_(W)  =  max 

^  w  •  k  =  1 

W2k»  k  ±, 


...  K(WV  'V  '  XT2Jk] 

C  +  hl(w“w2k) 


where 


f2^u2Jk’  u2j]2  111  PHjk  “  FH,"2jkI  U2jk^ 

T2jk  =  *  *2’  (w2Jk> 


3(W)  =  .  3Tl,  ....  s  [W  *  WV  *  h2(W-%>| 


The  dynamic  programming  methods  of  Chapter  II  of  reference  5  are 
next  applied  to  the  recurrence  relations  h^,  h^  and  h^  for  some  fixed 
X  to  obtain  a  solution  for  the  optimum  value  of  h^(w)  “  F*  and  for  the 
corresponding  optimum  values  of  w^8 9  and 

Next  the  corresponding  optimum  values  of  v^  and  v^  are  com- 

puted  using  the  computer  determined  optimum  values  of  9  and  < 

The  sum  (v^  +  v^  4*  v^)  is  then  compared  with  the  available  payload 
volume  V.  Probably  v^  +  v^  +  v^  f  V. 
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Then  if 


v!  +  v2  +  v3>  V 

the  value  of  X  must  be  increased  and  the  dynamic  program  repeated 
until  the  minimum  value  of  X  is  found  which  results  in 

V1  +  v2  +  v3  <V 
If  v^  +  v2  +  v^  <V 

the  value  of  X  is  decreased  until  the  minimum  value  of  X  is  found 
which  results  in 

V1  +  v2  +  v3  <V 

The  optimum  value  of  h^w),  w^,  w^'  and  w3  determined  for 

such  a  X  specify  the  solution  values  of  the  three  study  variables 0 
The  Lagrange-multiplier,  X  $  may  be  thought  of  as  a  cost 
coefficient .  If:  X  =  ° 

vl  +  v2  +  V3  ~V 


then  the  volume  restraint  may  be  thought  of  as  having  cost  nothing 0 
That  is,  the  measure  of  effectiveness  has  been  optimized  solely  with 
respect  to  the  weight  restraint  and  the  result  has  met  the  volume 
requirements*  On  the  other  hand,  if  X  -  o 

’l*"!*  v3>  v 

then  some  positive  value  must  be  assigned  to  X  which  will  result 
in  meeting  the  problem  volume  restraint  Vj^  +  +  v^^V. 

Ultimately  the  subsequent  process  of  optimization  with  some 
positive  value  of  X  will  result  in  smaller  optimum  values  of  one 
or  more  of  the  weight  variables  with  a  correspondingly  reduced  value 
of  the  measure  of  effectiveness. 
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This  reduction  in  the  solution  value  of  the  measure  of  effective- 


ness  may  be  considered  as  the  cost  of  having  a  volume  restraint 0  Of 
course  if  weights  and  volumes  are  interchanged  in  the  model  the 
resulting  X  will  be  a  measure  of  the  cost  of  the  weight  restraint 0 
If  the  solution  value  of  probability  of  kill  is  desired  then 


h^(w)  +  X  (vL  +  v2  +  v^)J 
,  v^,  v^  and  v^  are  solution  values. 


Pk  =  exp 


where  h^(w) 


3 


The  computer  flow  chart  to  accomplish  this  dynamic  programming 
procedure  is  contained  in  Appendix  B. 

It  may  be  of  interest  to  allow  the  values  of  the  available 
payload  weight  and  volume  (w  and  V )  to  vary  in  the  vicinity  of  their 
preset  values.  Then  the  sensitivity  of  the  effectiveness  measure ^ 
to  such  small  variations  may  be  examined  in  order  to  determine 
the  possible  value  of  modifications  to  the  weapon  delivery  systems  to 
allow  for  a  different  (more  effective)  load  capability „ 

Input  Requirements  for  the  allocation  model 
W  =>  Available  payload  weight  (lb) 

V  =  Available  payload  volume  (ft^) 


W  ,  W  =  lower  and  upper  limits  respectively  on  the  allowable  range 
of  fuel  weight  (lb) 

t  i  t 

W  ,  W  =  lower  and  upper  limits  respectively  on  the  allowable 


range  of  weight  for  modifications  to  the  propulsion 
system  to  provide  a  maximum  (attack)  speed  capability  in 
excess  of  a  minimum  preset  value  (lb) 
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I  I  I  I  I 

W  ,  W  =  lower  and  upper  limits  respectively  on  the  allowable 

range  of  explosive  weight  (lb) 

t  v 

gl(wl);  i  =  1,  2  ,  3  =  the  relationships  between  the  i  th  weight 

variable  and  the  i  th  volume  variable; 
i  v 

i  =  1,  2  ,  3 

Pa(Wi)  =  probability  of  target  acquisition  expressed  as  a  function 
of  the  fuel  weight  for  acquisition 

P„(w* )  w*  1  )  =  probability  of  hitting  the  target  given  that  target 

n  2  2 

acquisition  has  occurred  expressed  as  a  function  of  the 
fuel  weight  for  hit  and  of  the  weight  for  modifications 
to  the  propulsion  system  to  provide  a  maximum  (attack) 
speed  capability  in  excess  of  Sq 

P  (w  )  a  probability  of  target  kill  given  that  a  hit  has  occurred 

E  J 

expressed  as  a  function  of  the  explosive  weight 
Sq  =  minimum  preset  value  of  attack  speed  (yd/sec) 

£w  =s  size  of  the  increments  of  weight  of  the  weight  variables^ 
all  of  which  are  varied  in  increments  of  equal  size  (lb) 

The  payload  weight  and  volume  are  determined  by  the  capabilities 
of  the  weapon  delivery  systems.  The  limits  of  the  weight  variables 
are  determined  by  a  preliminary  study  which  reflects  limiting  design 
considerations  and  reduces  the  possible  ranges  of  the  weight 
variables  to  reduce  computer  time.  The  weight~volume  relationships 
are  determined  by  engineering  considerations  reflecting  the  state 
of  the  art.  The  minimum  attack  speed  is  determined  by  intelligence 
estimates  or  CNO  requirements.  The  size  of  the  weight  increments 
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is  determined  by  considerations  of  required  accuracy  of  results  and 
available  computer  time.  The  probabilities  of  acquisition,  hit  and 
kill  may  be  obtained  in  a  number  of  ways  including 

1.  analytical  derivations  in  terms  of  obtainable  or  estimable  inputs 

2.  simulation  (war  gaming) 

3.  extrapolation  of  empirical  data  from  operational  or  development 
tests 

4*  a  combination  of  the  above 

The  remaining  chapters  provide  samples  of  possible  analytical  methods 
of  providing  the  three  required  probability  inputs  which  are  specifi¬ 
cally  designed  to  be  compatible  with  the  allocation  model  requirements 
It  is  possible  to  extend  the  model  to  provide  for  cost  restraints 
The  procedure  would  require  a  second  Lagrange-multiplier  in  the 
function  to  be  maximized  (f1  )<>  Additional  inputs  required  would  be 
the  relationsuips  between  the  costs  and  the  corresponding  weight 
variables.  It  is  considered,  however,  that  cost  restraints  are  not 
significant  to  the  problem  considered  in  this  paper.  Variations  in 
the  amount  of  fuel  and  amount  of  explosive  have  little  effect  on  the 
costs  involved.  Even  changes  in  attack  speed  within  the  limits 
appropriate  for  consideration  will  probably  have  negligible  cost 
effect  when  compared  to  the  total  cost  of  developing,  producing  and 
operating  a  weapon  in  the  course  of  its  life  cycle. 
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CHAPTER  III 


PROBABILITY  OF  ACQUISITION 

Consider  a  tactical  situation  in  which  a  weapon  is  delivered  by 
some  air  delivery  weapons  system  to  a  splash  point  on  the  surface  of 
the  water  in  the  vicinity  of  a  submarine  target  and  make  the  following 
assumptions . 

Assumption  1;  The  splash  point  defines  the  center  of  the  top  of  a 
right  circular  cylinder  called  the  target  volume  cylinder.  The 
probability  is  quite  small  and  is  assumed  zero  that  the  target  is  out= 
side  of  the  target  volume  cylinder.  The  target  is  placed  at  random 
within  the  target  volume  cylinder. 

Assumption  2:  The  weapon  dives  to  its  assigned  search  depth  on  a 
straight  line  track  and  commences  a  right  circular  search  turn  at 
constant  depth,  speed,  and  turn  rate  immediately  upon  reaching 
search  depth.  The  weapon  searches  in  the  active  acoustic  mode 
commencing  its  first  ping  upon  commencement  of  the  search  turn. 

The  weapon  continues  its  turn  until  target  acquisition  or  fuel  ex= 
haustion. 

Assumption  3:  The  weapon  does  not  acquire  a  false  target. 

Assumption  k:  The  search  turn  of  the  weapon  is  approximated  by 

rotating  the  weapon  about  a  vertical  axis  through  the  point  at  which 
the  weapon  commences  search. 

Assumption  5:  While  the  weapon  is  diving  to  search  depth,  the  radius 

of  the  target  volume  cylinder  increases  at  a  rate  equal  to  the  target 

cruise  speed. 
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Assumption  6:  The  target  is  alerted  on  commencement  of  the  weapon  search 
turn  by  the  first  'ping'  of  the  weapon.  Thereafter  the  target  opens  at 
constant  depth  and  speed  radially  from  the  projection  of  weapon  position 
on  the  horizontal  plane  through  the  target a 

Assumption  7;  During  search  the  radius  of  the  target  volume  cylinder 
increases  at  a  constant  rate  equal  to  the  target  speed. 

Assumption  8:  The  'definite  range*  law  as  developed  in  Chapter  II  of 
reference  6  applies  to  acquisition.  Thus  on  a  given  radial  during  a 
ping  cycle  under  specified  conditions  there  is  a  range  beyond  which 
acquisition  probability  is  zero  and  within  which  acquisition  probability 
is  one.  Regions  closer  than  the  'definite  range'  limits  are  considered 

scanned  „ 

Assumption  9:  During  search  the  weapov  makes  a  series  of  independent 
glimpses  for  the  target,  one  for  each  ping  cycle. 

Assumption  10:  Target  cylinder  volume  is  considered  constant  during 
each  ping  cycle  at  the  value  valid  at  the  initiation  of  the  ping. 
Assumption  11:  The  water  medium  is  considered  to  be  isotropic. 

That  is,  all  environmental  parameters  that  affect  acoustic  performance 
have  a  constant  value  throughout  the  considered  volume. 

Assumption  13:  Consider  planes  containing  the  transducer  and  angles 
measured  at  the  transducer.  Consider  the  vertical  plane  containing  the 
weapon  longitudinal  axis  as  the  vertical  reference  plane.  As  the 
weapon  sweeps  through  an  angle  <P  in  the  horizontal  plane  it  sweeps 
through  a  lesser  angle  <p »  in  a  plane  perpendicular  to  the  vertical 
reference  plane  but  inclined  to  the  horizontal.  Also  if  the  inclination 
of  the  plane  inclined  to  the  horizontal  as  measured  in  the  vertical 
reference  plane  is  «  then  the  inclination  of  the  plane  inclined  to  the 
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horizontal  as  measured  in  a  vertical  plane  other  than  the  reference 
plane  is  a  lesser  angle  « ' .  In  this  paper  it  is  assumed  that 

©C  *  =  ©C 

In  practical  situations  likely  to  be  considered,  this  assumption 
results  in  much  smaller  errors  than  those  created  by  treating  the 
range  problem  in  discrete  steps  as  is  done  in  this  chapter. 

Define : 

P  =  p  (i)  =  P  (”  acquisition  in  I  ping  cycles  ] 

A  A  b 

where  P  f  xj  =  probability  of  event  X  happening 

g.  =  P  [" acquisition  during  the  i  th  ping  cycle]  j  i  -  1,  •••,  I 

w  =  weapon  fuel  expenditure  rate  at  search  depth  and  speed  (Ib/sec) 

Prr  =  weapon  active  sonar  pulse  repetition  rate  for  search 
s 

(pings/sec) 

w  =  fuel  used  by  weapon  to  reach  search  depth,  (lb) 

OSD 

w  =  fuel  available  for  acquisition,  (lb) 

Then,  as  developed  in  Chapter  2  of  reference  6, 

(1-gj) 


PA=1 
and  since 


i=l 


(W1  “  wOSD')  PRRs 


v 


It 


=  I 


and  X~  =  greatest  integer  —X  and  where  I 

is  the  number  of  complete  ping  cycles  available  for  the  event  of  acquisi^ 
tion,  it  may  be  said  that  P  is  a  function  of  wL  and  a  table  may  be 
computed  of  vs  wx  where  wL  is  allowed  to  vary  between  W  and  W  . 
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Of  course  if  I  is  less  than  1  then  PA  is  defined  to  be  zero*  The 
probability  of  acquisition  in  a  specified  glimpse  (g^)  con^ 

sidered  as  a  coverage  factor*  The  volume  containing  the  target  with 
positive  probability  is  a  region  to  be  covered c  A  portion  of  the 
region  to  be  covered  is  scanned  (covered)  by  the  weapon  in  accordance 
with  the  definite  range  law*  Because  the  target  has  been  randomly 
placed  in  the  target  volume  cylinder  (region  of  target  uncertainty 
to  be  covered)  then 

volume  scanned  on  i  th  ping  cycle 
®i  ~  target  volume  during  i  the  ping  cycle 

Thus  the  problem  is  reduced  to  finding  what  portion  of  the  volume 
of  target  uncertainty  is  scanned  by  the  weapon  on  any  ping*  This  will 
be  different  for  every  ping. 

The  following  development  for  is  divided  into  eight  consecutive 
steps *  The  development  is  amplified  as  appropriate  and  complete  flow 
charts  for  a  computer  program  are  provided  in  appendices  C  1  through  C  8* 
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1.  Available  Echo  Excess  (B^k) 

Consider  a  vertical  plane  containing  the  longitudinal  axis  of  the 
weapon  at  the  time  of  ping  transmission  Consider ,  as  in  Figure  1,  a 


second  plane  perpendicular  to  the  first  and  containing  the  weapon 
transducer  such  that  it  makes  an  angle  j  with  the  weapon  longitudinal 
axis,  A  set  of  such  planes  is  considered  such  that  the  angle  j  is 
associated  with  the  j  th  plane.  In  the  vertical  plane  angles  are 
measured  positive  counterclockwise  from  the  reference  axis. 

Design  or  operational  considerations  will  make  it  necessary  to 
prohibit  reception  of  signals  from  off  the  transducer  axis  in  the 
vertical  plane  in  excess  of  some  limting  angle  J,  not  greater  than 
90  degrees.  Allowing  one  degree  between  each  of  the  considered  planes  j^, 

the  index  j  will  assume  the  range  of  values  -J,  =J  +  1#  °  °  °i> 

J.  Thus  2J+1  planes  are  to  be  considered,  one  containing  the  trans~ 

ducer  axis  and  J  planes  on  either  side* 


Consider  increments  of  target  range  of  100  yards  in  the  plane  j. 

For  an  assumed  target  positioned  in  the  j  th  plane  at  a  range  of  100k 
yards  on  the  transducer  azimuth  axis  projection  in  the  plane  j  (point  P 
in  Figure  l),  there  exists  an  echo  excess  (Bjk)  by  which  the  signal  level 


<SV 


of  the  echo  reflected  from  the  target  exceeds  the  background 


masking  level  at  the  weapon  transducer  (M(BG)k).  That  iss 


>  o  »  K 
S> 


B^k  =  SL^k“  M(BG)k;  j=  -J,  •**>  J5 

Assuming  a  definite  range  law  of  acquisition,  a  negative  echo  excess 
yields  zero  probability  of  acquisition  and  a  positive  echo  excess  yields 
a  probability  of  acquisition  of  one.  The  maximum  range  of  acquisition  on 
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X  direction.  Line  ox  is  the  transducer  axis  extension. 


ox,  oP,  and  oy  are  in  the  vertical  plane.  Lines  oP  and  oz  are 
in  the  plane  j  .  The  polar  coordinates  of  the  point  P  in  the 
vertical  plane  using  the  x  axis  as  a  reference  are  (lOOkp  j). 
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the  azimuth  axis  in  the  j  th  plane  is  thus  the  range  of  1  yards 

which  B  equals  zero.  At  ranges  less  than  100K.  yards  there  is  an 
jk  J 

echo  excess  of  B above  that  required  for  acquisition  on  the  azimuth 

jk 

axis.  This  excess  is  available  for  use  off  the  aximuth  axis  d 

The  signal  level  varies  with  assumed  target  range  and  with  the 
elevation  angle  of  the  target  above  or  below  the  tranducer  axis  in 
accordance  with 

SLJk  =  WS  +  TS  -  2  Hfc  +  10  log  ^(J)  b2^] 

where : 

WS  =  on- axis  source  level  (db  at  1  yard  from  transducer 
TS  =  target  strength  of  target  (db) 

H  =  one  way  transmission  loss  at  a  range  of  100k  yards  (db) 
k 

t»i(j)  =  normalized  vertical  pattern  function  for  transmitting 
pattern  (i.e.,  the  ratio  of  intensity  at  the  angle  j 
off  the  vertical  axis  but  on  the  azimuth  axis  to  intensity 
on-axis  in  both  the  vertical  and  aximuth  planes, 
b  (j)  =  normalized  vertical  pattern  function  for  receiving  pattern. 
There  are  several  different  kinds  of  masking  noises  and  it  is  their 
total  effect  which  is  considered  as  the  interference  background  masking 
level.  Included  are  self-noise,  target  noise,  ambient  noise,  and 

volume  and  boundary  reverberation. 

There  is  also  an  enabling  level  designed  into  the  weapon  electrical 
circuitry  to  protect  the  acoustic  system  from  excessive  interference 
and  to  avoid  interpretation  of  the  interference  as  a  signal.  As  used 
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here,  the  weapon  enabling  level  (WEL^)  includes  the  weapon  recognition 
differential.  It  is  a  function  of  range  (or  time)  and  decreases  with 
increasing  range.  The  background  masking  level  (m(BG )r )  is  the  larger 
of  the  two:  interference  background  masking  level  and  the  weapon 
enabling  level. 

The  output  of  this  step  is  a  table  of  the  available  echo  excess 

(BJk;  i=  -J,  •••,  J;  k-  l-  •••>  K')>- 
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2. 


Scanned  Volume  Radials;  General  Case 
In  each  of  the  planes  j  there  is  now  defined  for  each  value  of 
range  (100k)  an  available  echo  excess  B^k  for  acquisition  on  the 
azimuth  axis.  This  available  excess  may  be  used  to  find  limits  for 
off-axis  acquisition.  Define  an  angled  to  be  the  angle  through 
which  the  weapon  turns  while  a  ping  travels  to  a  range  of  100k  and 


returns , 
Then 


*k  = 


200 ka)g 


where 

0)  =  weapon  search  turn  rate  (deg/sec) 

c  =  velocity  of  sound  in  the  medium  (ft/sec) 

For  echo  ranging  In  the  plane  j  to  a  range  o£  100k»  an  energy  pulse 

transmitted  at  angle  p  will  return  at  the  angle  p  -*k- 

The  echo  excess  Ajfc  ( p)  required  to  acquire  the  target  at  range  100k 

in  plane  j  at  angle  p  off  the  aximuth  axis  is 

A  k(p)  =  -  10  log  b3  (p)  -  10  logb^CP  k) 
where  b3  (  p)  =  normalized  aximuth  pattern  function  for  transmitting 

pattern 

(  p)  =  normalized  aximuth  pattern  function  for  receiving 
pattern 

If  the  required  echo  excess  A Jfe  ( f»  exceeds  the  available  echo  excess 
B  k  then  acquisition  will  not  take  place  at  range  100k  at  an  azimuth 
removed  p  degrees  from  the  reference  axis  of  plane  j.  However  if 
A  (  p)  does  not  exceed  Bjk,  acquisition  takes  place. 
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jn  each  plane  j  and  at  each  target  range  100k  there  are  two  limit ing 
angles  between  which  acquisition  will  occur  with  probability  one  and 
outside  of  which  there  will  be  no  acquisition.  Looking  from  the  trans¬ 
ducer  toward  the  potential  target  define  the  set  of  limiting  angles 

8  with  m  1  for  the  limiting  angle  to  port  and  m  =  2  for  the 
Kjmk 

limiting  angle  to  starboard.  As  shown  in  Figure  2f  the  reference 
axis  for  polar  coordinates  in  the  j  th  plane  is  the  projection  of 
the  transducer  axis  at  the  time  of  ping  transmission  in  that  plane. 
Positive  angles  are  measured  clockwise  from  the  reference  axis  in  all 
planes  j  and  in  any  horizontal  plane.  The  pattern  may  not  be  centered 
on  the  axis  and  thusp^  may  be  positive  or  negative.  Design  or 
operational  considerations  will  make  it  necessary  to  prohibit  reception 
of  signals  from  off  the  transducer  axis  in  the  azimuth  plane  in  excess 
of  some  limiting  angle  a  in  the  direction  opposite  to  the  search  turn 
and  «  '  in  the  direction  of  the  search  turn.  This  restricts  the 
angles  P  as  follows 

Pjlk  “  ^k 

Pj2k  “  ^  k  |  “  * 

Recall  that  in  plane  j  the  maximum  range  for  a  non- negative  Bjk 
was  computed  to  be  100K^  yards.  For  a  straight  running  weapon  this 
would  be  the  maximum  acquisition  range  in  plane  j.  Actually  an 
acquisition  range  of  100K'  yards  will  not  be  realized  in  plane  j  in 
this  problem  because  the  weapon  is  turning  at  some  rate.  The  energy 
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100k 


O  o 
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Figure  2:  Acquisition  limiting  angles  in  the  plane  j  „  Line  OP  is  the  reference  axis 
The  x  axis  is  the  transducer  axis  projection,,  In  this  case  j2k  is  zero 


from  a  given  source  direction  will  never  return  as  an  echo  in  the  same 
direction.  Thus  there  will  always  be  some  penalty  for  off~axis  trans- 
mission  or  reception. 

Unless  the  turn  rate  is  extremely  high>  the  insonified  volume  of 
adjacent  weapon  pings  will  overlap.  If  a  given  ping  scans  part  of  the 
volume  insonified  by  the  previous  ping>  the  overlap  cannot  be  credited 
to  the  volume  scanned  by  the  later  ping  because  the  overlap  volume  was 
removed  from  the  volume  of  target  uncertainty  by  the  previous  ping* 

The  concept  of  increased  probability  of  detection  because  of  overlap 
is  not  valid  in  this  case  because  of  the  use  of  the  definite  range 
law.  Thus  all  overlap  volume  should  be  credited  to  the  ping  which 
first  scanned  it.  This  would  give  all  pings  subsequent  to  the  first 
ping  during  the  first  search  turn  an  equal  scanned  volume  less  than 
that  of  the  first  ping.  In  this  paper  the  situation  is  approximated 
by  dividing  the  volume  of  overlap  between  the  pings  concerned  to  give 
all  first  turn  pings  an  equal  volume  of  insonif ication ,  This  approxima° 
tion  decreases  P  (l)  but  the  error  in  P  (i)  decreases  with  increasing 

n  n 

I  until  the  end  of  the  first  search  turn  when  it  becomes  negligible. 

The  plane  of  maximum  overlap  is  plane  0.  In  this  plane^  for  small 
k,  f*Q2k  for  PinS  i  will  be  larger  than  b  +  p  for  PinS  (i+l). 

See  Figure  3*  Here  b  is  the  angle  between  the  axes  of  pings  i  and 
(i+l)  or  the  angle  through  which  the  weapon  turns  between  pings.  As  k 
becomes  larger  the  differences  between  P  ^  for  ping  i  andb  +  p^^  f°r 
ping  (i+l)  become  smaller.  Define  pf  to  be  the  angle  at  which  the 
two  are  equal.  Thus  the  right-most  angle  for  acquisition  in  any  plane 


3^ 


Figure  3:  Overlap  of  insonified  volume  between  adjacent  pings. 
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is  not  greater  than  p'  and  the  leftmost  angle  is  not  less  than  ?'-*>• 
Some  of  the  extreme  planes  j  may  not  be  useful  for  acquisition  at  all. 
By  including  for  consideration  only  those  planes  in  which  acquisition 
may  take  place  the  index  j  takes  on  the  new  range  from  J1  below  the 
axis  to  J1 1  above  the  axis. 

The  range  or  radial  distance  associated  with  the  angle  pjmk  is 
100k  yards.  See  Figure  4.  There  is  a  range  maximum  called  the  ping 
interval  range  defined  by  the  time  between  pings  (l/PRR,)  and  the 
velocity  of  sound  in  the  medium  (c).  Further,  for  planes  inclined 
above  the  horizontal  there  is  another  range  maximum  defined  by  the 
range  at  which  the  radial  at  angle  Pjmk  breaks  the  surface  of  the  water 
Similarly  for  planes  inclined  below  the  horizontal  there  is  a  range 


maximum  defined  by  the  range  at  which  the  radial  at  angle  pjmk  breaks 
the  maximum  target  depth  (D^).  Volume  scanned  below  the  maximum 
target  depth  does  not  contribute  to  acquisition  probability.  Thus  the 
acquisition  range  at  any  angle  Pjmk  is  the  minimum  of  the  ranges  just 
discussed  to  be  denoted 

Since  the  index  k  is  no  longer  of  value  to  indicate  range  it  is 

useful  now  to  relabel  the  radials  such  that  the  index  k  indicates  the 

separation  from  the  leftmost  radial  in  the  plane.  See  Figure  5.  The 

leftmost  radial  is  labeled  (R^,  Pjl).  The  general  radial  is 

(R  B,  )•  The  rightmost  radial  is  (R.r  ,  )•  For  Planes  close 

'  jk7  rjk  J  j  J  j 

to  the  axis,  probably 


PjK.  -  p  '  and  pjK_  ~  Pjj.  -  b 
This  may  not  be  true  of  the  extreme  planes. 
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Figure  4:  Acquisition  range  at  angle  Line  OR  has  length 

100k.  Line  SS1  denotes  the  surface.  Length  01  is  the 
ping  interval  range.  Lined  area  is  scanned. 


37 


Figure  5:  Relabeled  radials  in  plane  j„ 
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The  output  of  this  step  is  tables  of 

f>jk and  V  3  =  J’’  J":  k  = 

plus  the  angles  p  and  b  • 


K. 

J 
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3.  Scanned  Volume:  General  Case 

The  radials  computed  in  step  2  define  areas  of  insonif icat ion  in 
each  of  a  set  of  planes  j.  The  area  in  each  plane  may  be  computed  as 
the  sum  of  a  set  of  triangles  as  shown  in  Figure  6. 

y1 

Aj  =  (1/2)  pp  Rjk  Rj(k+1)  Sin  (pj(k+l)  “  Pjk} 


The  irregular  shape  of  the  areas  in  the  planes  j  makes  it  difficult 
to  work  with  them  in  their  present  form.  Define  a  set  of  isosceles 
triangles  (one  for  each  plane  j)  such  that  the  triangle  areas  are 
equal  to  the  scanned  areas  A  in  planes  j  and  the  common  apex  angle 


is  the  angle  turned  by  the  weapon  per  ping  (b). 
altitude  of  each  triangle  is 

f  "l1/2 


H 


3 


tan(b/2) 


See  Figure  7.  The 


The  set  of  triangles  so  defined  approximates  the  areas  scanned  and 
makes  it  possible  to  obtain  the  total  scanned  volume  SVQ1  on  ping  zero 
on  the  first  turn.  Ping  zero  is  a  ficticious  ping  denoting  the 
general  case.  The  computation  of  SV^  (scanned  volume  on  i  th  ping  of 
first  search  turn)  is  identical  to  that  for  SVQl  to  be  developed  next. 

Define:  p  =  weapon  angle  of  attack  (pitch)  (deg). 

The  plane  j  makes  an  angle  of  (j+p)  degrees  with  the  horizontal. 

The  set  of  isosceles  triangles  is  used  to  define  a  set  of  pyramids. 
A  vertical  cross-section  of  the  set  of  pyramids  is  shown  in  Figure  8 
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Figure  6:  Area  of  insonif ication  in  plane  j „ 


hi 


Figure  7: 


Equivalence  of  area  and  isosceles  triangle 
OHH1  of  altitude  and  apex  angle  b  . 
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ping  zero,  x  is  the  horizontal  axis,  y  is  the 


vertical  axis,  and  OQ  is  the  transducer  axis 
extension. 


^3 


and  a  single  pyramid  is  shown  in  Figure  9* 
pyramid  is 

V.  =  1/3  l”area  of  base  altitude 


The  volume  of  a  single 


As  shown  in  Figure  9>  2X  is  the  average  value  of  one  of  the 
base  dimensions  and  Y  is  the  value  of  the  other  dimension. 

Thus  the  base  area  is  2XY. 


X  =  |~(Hj  +  tan(b/2)]/2 

Since  sintf>  =  h/faj  =  sin  1  )/l 


the  altitude  is 

h  =  (Hj  sin  1°  )A 

and 

Vj  =  (l/3)(2XY)h  =  (l/3)(Hj+H(j+1))  H(j+l)  sin  1®  tan(&/2) 


The  volume  scanned  is 


J"-l 


SVrt,  = 


’01 


j=J' 


The  output  of  this  step  is  tables  of  H^j  j  =  J* 9  •••,  J* 1 
and  the  volume  SVq^. 


kk 


1*5 


Figure  9s  Pyramid  volume  V  .  c  Line  OE  is  0  Line  OH11  is  H( 


4.  Special  Volume  Elements:  General  Case 

For  the  case  of  search  at  constant  depth  in  a  continuous  circle 
as  is  being  considered  in  this  paper,  the  situation  changes  drastically 
upon  completion  of  the  first  search  turn  circle 0  Subsequent  search 
turns  will  effectively  scan  much  less  volume  because  the  weapon  retraces 
its  path  of  the  previous  turn.  Define  the  maximum  available  number  of 
search  turns  (U)  to  be 


U  = 


WOSD 


)  PRR  /M 
s 


+  1 


where  M  is  the  number  of  ping  cycles  in  one  search  turn 


M  =  ||36o/&||  +1 


IH' 


=  greatest  integer  _  x 


x||  ==  greatest  integer  <  x 

_ 1 

W  =  maximum  allowable  weapon  fuel  (lb) 

w^t  =  fuel  expenditure  rate  at  search  depth  and  speed  (lb/sec) 
Wqsd  =  fuel  used  in  dive  to  search  depth  (lb) 


(i/3)  Dyg  vQt 


v( 


OSD  sln^g 


l/3  Duc  =  weapon  search  depth  (yd) 

W  b 

=  average  weapon  fuel  expenditure  rate  in  dive  to  search 


<P 


DS 


depth  (lb/sec) 

=  weapon  speed  in  dive  to  search  depth  (yd/sec) 
=  angle  of  dive  to  search  depth  (deg) 


k6 


Consider  a  side  view  or  vertical  cross-section  of  the  scanned  volume 
of  a  single  ping  as  in  Figure  10c  The  vertically  lined  area  is  that 
scanned  by  the  weapon  on  the  i  th  ping  of  the  first  search  turn* 

Recall  assumption  6  that  the  target  opens  radially  from  the  weapon 
position  and  assumption  12  about  weapon  orientations  on  different 
search  turns.  Define  a  distance  of  TES  yards  as  the  distance  traveled 
by  the  target  during  the  time  it  takes  the  weapon  to  complete  one 
search  turn.  The  weapon  scanned  area  on  the  i  th  ping  of  the  second 
search  turn  is  equivalent  to  that  of  the  areas  lined  sloping  downward 
to  the  right  on  Figure  10.  The  weapon  actually  scans  the  same 
geographic  volume  on  each  search  turn  but  since  the  target  is  assumed 
to  open  radially  the  situation  is  equivalent  to  translating  the 
area  scanned  on  the  first  turn  a  distance  TES  to  the  left  for  each 
subsequent  turn.  If  the  target  is  in  the  area  lined  sloping  downward 
to  the  left  on  the  i  th  ping  of  the  first  search  turnp  he  will  escape 
prior  to  the  i  th  ping  of  the  second  search  turn.  Similarly  for  sub¬ 
sequent  search  turns.  A  target  above  the  critical  depth  CDA  or  below 
the  critical  depth  CDB  may  escape  but  a  target  at  any  intermediate  depth 
will  be  acquired  eventually  if  fuel  supply  permits  because  of  the  assump¬ 
tion  of  an  opening  target.  The  search  turns  after  the  first  are  desirable 
because  of  the  possibility  of  acquiring  a  target  1  trapped0  inside  the 
insonified  volume  of  the  first  search  turn.  Treatment  of  scanned 
volumes  on  the  second  and  subsequent  turns  will  differ  depending  on 
whether  the  volume  is  above  or  below  the  depths  CDA  and  CDB.  Although 
the  areas  scanned  above  depth  CDA  on  search  turns  1  and  2  are  equal  in 
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Figure  10;  Vertical  cross- section  of  a  single  ping.  Lines  xx“  have  length  TES 
Lines  yyc  and  lines  xx#  have  length  L  . 


Figure  10,  the  volumes  represented  by  these  areas  are  not  equal*  A  top 
view  of  ping  i  is  a  pie  shaped  slice  which  makes  the  volume  associated 
with  the  first  search  turn  area  larger  than  that  of  the  corresponding 
second  search  turn  area.  Thus  each  volume  element  must  be  calculated 
independently  and  the  volume  scanned  by  the  weapon  on  ping  i  of  search 
turn  u  (SV^)  is  the  sum  of  several  independently  calculated  volumes* 

In  this  step  the  volume  elements  in  the  general  case  (i-0)  are  computed* 
These  are  modified  as  necessary  for  each  ping  orientation  in  a  later 
step.  Although  Figures  11  and  12  and  the  development  to  follow  involve 
the  volumes  above  the  search  depth,  the  development  for  volumes  below  the 
search  depth  is  identical. 

Consider  the  vertical  cross-section  of  Figure  10*  The  cross- 

section  has  been  split  into  several  regions  or  volume  layers  identified 

on  the  left  margin  by  their  bottom  boundary  labeled  such  that  the 

bottom  boundary  of  a  level  (or  layer)  at  or  above  depth  CDA  has  a 

length  of  L,  inside  the  first  turn  scanned  region*  The  planes  j  are 
a 

considered  using  the  plane  (j  =  TJ)  of  maximum  H.  |sin  ( j  +  p)  |  as  a 
reference  plane.  In  Figure  10  below  search  depth 
HXJ  | sin  (TJ  +  p)|  =  H(TJ+1)|sin  (TJ+l+p)  |  . 


In  this  case  where  more  than  one  of  the  planes  j  has  the  same  maximum 

|  sin  (j  +  p)|,  the  reference  plane  is  chosen  to  be  that  with  the 

least  Hj  cos  (j  +  p).  Above  search  depth,  each  plane  j  such  that 

H.  sin  (i  +  p)<D„„  -  CDA  and  i>-TL  defines  a  level  L , .  Also  each 

j  J  WS  J  d 

plane  j  such  that  sin  ( j+  [1)^0^=  CDA  defines  a  level  L^„ 
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Referring  to  Figure  11,  the  volume  elements  LADV^  ,  MADV^,  and 
RADV^  (d  -  1,  d;  u  -  2,  eee,  U)  are  significant.  The  subscript  d 

indicates  the  level  and  the  subscript  u  indicates  the  search  turn 
number.  The  level  L-  is  the  final  level  above  the  critical  depth  CDA 
and  =  TES .  The  corresponding  volumes  below  the  search  depth  are 


LBDVdu'  mDVdu>  and  RBDVdu  (d  =  1>  R'>  u  =  2, 


du 


du 


U  ).  The  scanned 


volume  in  level  above  the  search  depth  on  the  u  th  search  turn  is 

ASV  =  LADV ,  +  MADV,  +  RADV, 

du  du  du  du . 


One  or  more  of  these  volume  increments  may  be  zero  as  is  LADV^ 

in  Figure  11.  The  corresponding  scanned  volume  below  the  search  depth 

is  BSV,  . 
du 

Referring  to  Figure  12, 

the  volumes  AMV^;  d  =  D+l;  •••,  D*  j  u  =  1,  U  are  significant. 

The  corresponding  volumes  below  the  search  depth  are  denoted 
BBVdu  (d  =  D  +  1,  D';  u  =  1,  U  ) 

For  levels  L,  below  depth  CDA  and  above  search  depth  and  for  the  u  th 
d 

search  turn,  AMV^  is  the  trapped  volume  of  target  uncertainty  remain¬ 
ing  after  the  ping  cycle  is  complete.  The  volume  scanned  in  level 
above  the  search  depth  on  the  u  th  search  turn  is  thus 


ASV  =  AMV  ,/  .  N  -  AMV  , 

du  d^u-1;  du 

The  corresponding  scanned  volume  below  the  search  depth  is  BSV^. 

In  this  step  all  calculations  are  done  for  ping  zero  and  the 
quantities  computed  for  ping  zero  are  saved  to  be  modified  later  as 
appropriate  for  the  i  th  ping.  The  total  scanned  volume  for  ping 
zero  of  the  u  th  search  turn  is 
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(T-P) 
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(T-P) 
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D'  D® 

STou=  2  “vdu  ♦  jr  Kvdu 

d=l  d=l 

The  output  of  this  step  is  the  number  of  pings  (m)  in  a  search 
turn,  the  maximum  available  number  of  search  turns  U,  several  tables 
of  parameters  described  in  appendix  Ck  used  in  computing  the  special 
volumes,  and  tables  of  the  following  special  volumes  for  ping  zeros 


LADV ,  .  MADV,  ,  RADV,  :  d  =  1, 
du7  du7  du7  7 

u  =  2, 

°  j  u 

LBDVdu,  MBDVdu,  RBDV^J  d  .  1, 

•  *  ’} 

u  «  2, 

°  °  $  u 

AHVdu;  d  =  D+  1,  •••,  D'j  u  . 

lf 

U 

BMVdu;  d  -D+  1,  •••,  D'j  u  = 

1  •  •  0 

L  9  9 

U 

ASVdu;  d  .  1,  •••,  D’j  «  .  2, 

•  • 9  9  u 

BSVdu;  d  *  l>  •••»  5.' ;  “  -  2. 

•••,  u 

SV0u:  u  *  2<  U 
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5.  First  turn  scanned  volume:  i  th  ping 

Throughout  the  search  phase,  the  weapon  is  displaced  from  the 
center  of  the  target  volume  cylinder  a  horizontal  distance  R^ .  The  base 
radius  of  the  target  cylinder  is  some  distance  at  the  time  the 

weapon  splashes  into  the  water.  The  radius  then  increases  at  a  rate 
(Stc)  ec*ua^  to  t*ie  target  cruise  speed  to  a  distance  R7  at  the  time 
of  the  first  weapon  ping.  Thereafter  the  radius  increases  at  a  rate 

(S  c)  equal  to  the  alerted  target  escape  speed  and  has  a  value  R. 

1  b  i 

upon  commencement  of  the  i  th  ping  cycle. 

In  the  situation  of  Figure  13,  a  top  view  of  the  target  cylinder 
of  increasing  volume  with  the  first  few  ping  orientations  and  scanned 
volumes  in  the  first  search  turn,  the  insonified  volume  of  ping  1 
includes  a  large  volume  outside  the  target  cylinder.  Such  volume  does 
not  contribute  to  acquisition  probability  because  there  is  an  assumed 
zero  probability  that  the  target  is  outside.  Therefore  the  total 
insonified  volume  must  be  reduced  to  the  useful  scanned  volume  inside 
the  target  cylinder.  In  the  situation  depicted,  the  amount  of  dis~ 
carded  volume  is  reduced  with  each  ping  until  ping  9  when  all  of  the 
insonified  volume  is  within  the  target  volume  cylinder. 

The  projection  of  the  radial  (r^*  p^)  pl&ne  j  to  the  hori=> 

zontal  plane  through  the  transducer  is  the  line  OP  of  Figure  Ik 
denoted 

where 

R,jk -  Rjk  c°s  > 


5^ 


2 


Figure  13:  Top  view  of  target  cylinder  of  increasing 

volume  with  the  insonified  volumes  of  several 
sample  pings  of  the  first  search  turn  depicted « 
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Figure  14 :  The  projection  of  plane  j  in  the  horizontal  at 

search  depth.  Line  00'  is  the  plane  j  reference 


axis  projection.  Line  IP  is  CR^- 
r'.,  .  Line  IP'  is  R,  .  Line  OP*  is  R* 


Jk* 


Line  OP  is 

i 

ijk° 
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Use  of  the  law  of  cosines  enables  determination  of  the  distance 


I  P  denoted  CR^*  This  distance  is  compared  to  the  radius  (R^)  of  the 

target  cylinder  on  the  i  th  ping.  If  the  radius  (R^)  smaller^,  it 

defines  the  new  triangle  IOP1  and  the  projection  of  the  side  0Pff  back  to  the 

plane  j  defines  a  new  acquisition  range  (R  ,)  to  be  associated  with  the 

l  jk 

angle  p.,  .  That  is 

J  t 

Rijk 


R. 


ijk  cos 


ijTTi 


However  if  the  distance  CRt1  is  the  lesser  of  the  two  then 

jk 


R.  =  R., 
ijk  jk 


The  remaining  portion  of  the  computation  for  scanned  volume  (SV^^) 
on  the  i  th  ping  of  the  first  search  turn  is  identical  to  that  in 
step  3  for 

The  output  of  this  step  is  the  table  of  volumes  SV^j  1  =  1^  M=X 

and  the  distances  R  and  R.:  i  =  1,  °°°f  M-l. 

o  1 
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6.  Scanned  volumes  subsequent  to  first  turn*  i  th  ping 

Just  as  the  first  turn  scanned  volume  had  to  be  corrected  for  the  1  th 
ping  because  of  the  effect  of  the  target  cylinder  radius^  so  must  the 
scanned  volumes  on  subsequent  turns* 

Define  the  intersection  of  the  ping  i  bisector  with  the  target 
cylinder  radius  at  transmission  of  ping  i  to  be  a  distance  CH  from  the 
weapon . 

Referring  to  the  vertical  cross-section  of  the  insonified  volume 
in  Figure  15,  verticals  through  four  of  the  possible  intersections  have 
been  plotted  and  labeled  CHI  through  CH4 .  The  effect  of  each  of  these 
four  cases  will  be  discussed  as  it  applies  to  a  level  above  and  a  level 
below  the  critical  depth  CDA  (i,ej  levels  2  and  4)*  The  treatment  of 
the  situation  below  the  search  depth  is  identical* 

In  the  case  of  CHI  the  first  turn  scanned  volume  is  reduced* 

However  the  scanned  volumes  for  all  subsequent  turns  are  the  same  as 
those  computed  in  step  four  for  the  general  ping*  This  is  because 
CHI  lies  to  the  right  of  all  areas  representing  equivalent  volumes* 

In  the  case  of  CH2  there  is  still  no  reduction  for  volumes  in 
level  2*  However  there  is  a  small  reduction  in  the  volume  scanned  on 
the  second  search  turn  in  level  4,  The  volume  reduction  is  that 
volume  associated  with  the  area  in  level  4  to  the  right  of  CH3o 

It  is  computed  as  a  reduction  in  AMV^j  (the  volume  missed  in 
level  above  search  depth  on  the  first  search  turn)* 
ftuiS  (ASV22)plng  t  equals  (ASV22)plng  0  but 

(ASV>ping  i  is  Uss  than  (“Wng  O' 


58 


CH4  CHI  CH2  CHI 
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as  numbered o  Escape  volumes  have  primed  numbers  according  to  the  search  turn  of  escape 0 


In  the  case  of  CH3  there  is  a  reduction  in  the  scanned  volume  on 

the  second  search  turn  in  both  levels  and  0  In  level  the  volume 

(ASV^)  .  .  is  less  than  (ASV^)  .  A  by  an  amount  equal  to 

v  22  'ping  i  v  22  'ping  0  J  M 

RADV22  plus  a  portion  of  MADV^* 

In  level  the  reduction  still  amounts  to  a  portion  of  AMV^ 0 

In  the  case  of  CRk  there  is  no  scanned  volume  in  level  above 
search  depth  on  the  second  turn,  That  is  ASV ’  s  0  for  ping  i„ 

Further  there  is  a  small  reduction  in  level  in  the  scanned  volume 
for  the  third  search  turn.  For  level  there  remains  some  scanned 
volume  on  the  second  turn  and  there  is  also  a  reduction  in  scanned 
volume  for  the  third  turn. 

The  output  of  this  step  is  a  table  of  scanned  volumes 
SViu;  i  =  1,  •••,  M-lj  u  =  2,  U* 
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7-  Target  Volume:  i  th  ping 

Consider  that  portion  of  the  target  cylinder  volume  computed  at 
a  base  radius  of  within  the  radials  bounding  the  ith  ping*  See 
Figure  16*  The  area  contained  between  the  boundary  radials  of  the 
i  th  ping  and  within  the  cylinder  of  radius  multiplied  by  the 
height  of  the  cylinder  is  this  volume *  The  probability  that  the 
target  is  in  this  i  th  volume  slice  of  the  target  cylinder  is  the 
volume  of  the  slice  divided  by  the  total  cylinder  volume 0  This 
probability  (PQ^;  i  =  1,  •••,  M-l)  remains  constant  throughout  the 
problem  because  of  the  opening  target  assumption  number  six* 

From  Figure  16: 

pq  _  area  of  horizontally  lined  region 

1  7 T  Rj2 

The  target  volume  (TV^.;  i  =  1,  M-l)  associated  with  the  i  th 

ping  on  the  first  search  turn  lies  between  the  same  radial  limits  but 
within  a  target  cylinder  of  radius  instead  of  R^.  That  is 

TV^  =  (area  of  vertically  lined  region) 

Figure  16  shows  the  angle  (  b'  )  assigned  to  the  M  th  ping  to  be  less 
thanb  .  It  is 

b  '  =  360  -  (M-l)b 

This  smaller  angle  is  required  to  complete  the  search  turn  so 
that  the  i  th  ping  of  the  u  th  turn  is  oriented  in  the  same  direction 
as  the  i  th  ping  in  the  (u  +  1)  th  turn.  Thus  the  M  th  ping  is  in 
part  ficticious. 

The  output  of  this  step  is  thus  the  tables  of  TV^j  i  =  1,  •••,  M-l 
and  of  PQ^;  i  =  1,  •••,  M-l  and  the  angle  b*  . 
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Figure  16;  Top  view  of  target  cylinder  of  increasing  volume  . 

Lines  CT  and  Cl'  are  boundary  radials  of  ping  i. 

2 

Horizontally  lined  area  is  PQ^  times  TTR^ „ 
Vertically  lined  area  times  cylinder  height  is  TV^» 
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8.  Acquisition  Probabilities:  I  pings 

The  instantaneous  probability  of  acquisition  on  the  i  th  ping  is 
simply  the  ratio  of  the  volume  effectively  scanned  to  the  volume  of 
target  uncertainty. 

Define:  g^  =  P  ^acquisition  during  the  i  th  ping  cycle  of  the  u  th 

search  turn^j 

Then  glu  =  SVlu  ^  J  i  =  1,  •••,  M-i;  u  =  1,  •••,  U 

Recall  that  PQ^  is  the  probability  that  the  target  is  within  the 
radial  limits  of  ping  i. 

The  volume  TV^  is  the  total  volume  of  target  uncertainty  within 
the  radial  limits  of  ping  i,  The  volume  SV^  is  the  volume  scanned 
on  the  first  turn.  The  volumes  scanned  on  subsequent  turns  are  volumes 
measured  inside  the  first  turn  scanned  volume  and  under  the  conditions 
existing  for  the  first  turn.  These  volumes  are  then  moved  intact  with 
the  opening  taget .  Thus  the  applicability  to  all  turns  of  the  volume 
TV^  measured  using  the  first  turn  cylinder  radius. 

The  instantaneous  probability  of  detection  on  the  M  th  ping  is 
approximated  by 

«H>  =  *(M-l)u  {b,/b)  *  gl(u»l)  (&  *  b')/b 

The  M  th  ping  is  alloted  the  same  time  interval  as  any  other  ping0 
Errors  involved  should  be  small. 


63 


There  remains  the  requirement  to  relabel  the  to  be 

g^;  i  =  1,  •••,  I,  to  compute  I  as  a  function  of  as  described  in 
the  beginning  of  the  chapter,  to  compute  P  as  a  function  of  and 

A 

to  tabulate  the  values  of  versus  as  is  allowed  to  vary  from 

i 

W  to  W  .  This  final  table  is  the  input  to  the  dynamic  programming 
model . 
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9.  Input  Requirements 

A.  Medium 

c  =  velocity  of  sound  in  water  (yd/sec ) 
a  =  attenuation  coefficient  (db/yd) 

10  log  n>s  =  surface  reverberation  scattering  coefficient  (db) 
10  log  my  =  volume  reverberation  scattering  coefficient  (db) 

B.  Target 

TS  =  target  strength  (db) 

TN  =  radiated  noise  level  of  the  target  (db  at  1  yd) 

D_.  =  maximum  target  depth  (ft) 

S  =  target  speed  during  weapon  search  (yd/sec ) 

i.  j 

S  =  target  cruise  speed  (prior  to  becoming  alerted  on 
the  first  ping  cycle)  (yd/sec) 

C .  Weapon 

WS  =  transmitted  on-axis  source  level  of  weapon  (db  at  1  yd) 
AT  =  ping  length  (sec) 

PRR  =  weapon  active  sonar  pulse  repetition  rate  tor 
s 

search  (pings/sec) 

A  f  =  receiving  bandwidth  (cycles/sec J 

WEL^  =  weapon  enabling  level  tor  a  range  of  look  yds  (dbj 
SN  =  equivalent  isotropic  self°noise  level  at  search 
speed,  depth,  and  turn  rate  (dbj 
RD  =  recognition  differential  for  50$  probability  of 
acquisition  (db) 
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b ^  =  normalized  pattern  function  in  the  vertical  plane 

for  transmitting  pattern  (a  numeric) 

=  normalized  pattern  function  in  the  vertical  plane  for 
receiving  pattern  (a  numeric) 

b^  =s  normalized  pattern  function  in  the  azimuth  plane  for 
transmitting  pattern  (a  numeric) 
b^  =  normalized  pattern  function  in  the  azimuth  plane  for 
receiving  pattern  (a  numeric) 

J  =  limiting  angle  for  off=axis  reception  in  the  vertical 
plane  (deg) 

cr»  -  limiting  angles  for  off=axis  reception  in  the  azimuth 
plane  (deg) 

DR  =  directivity  index  for  receiving  (db) 
p  =  angle  of  attack  (pitch)  of  weapon  (deg) 

0)c  =  search  rate  of  turn  (deg/sec) 

^ DS  =  angle  of  dive  in  diving  to  search  depth  (deg) 

D  =  weapon  search  depth  (ft) 

W  j 

=  weapon  speed  in  dive  to  search  depth  (yd/sec) 

ST,„  =  weapon  speed  during  search 

WS 


W  =  maximum  allowable  weapon  fuel  weight  (lb) 
w  =  weapon  fuel  expenditure  rate  at  search  depth  and 
speed  (lb/sec) 

w =  average  weapon  fuel  expenditure  rate  in  dive  to 
search  depth  (lb/sec) 
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D.  Delivery  weapons  system  capability 

Rj  =  radius  of  target  volume  cylinder  at  time  of  weapon  splash 
Statistical  averages  and  distributions  are  available  for  the  medium 
inputs.  Intelligence  reports  and  estimates  or  CNO  requirements  may 
provide  the  target  inputs  except  for  the  target  speeds  which  may  be 
selected  in  the  higher  ranges  in  accordance  with  criterion  one 0 
The  weapon  inputs  are  available  from  weapon  designers  at  least  in  the 
form  of  educated  estimates.  Finally  the  input  may  be  estimated 
from  observations  of  the  performance  of  available  delivery  weapons 
systems . 
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CHAPTER  IV 


PROBABILITY  OF  HIT 

Consider  a  tactical  situation  in  which  a  weapon  acquires  a  submarine 
target  at  some  initial  slant  range  and  make  the  following  assumptions * 
Assumption  1:  The  weapon  makes  an  attack  and  a  series  of  independent 
reattacks  on  the  target  culminating  in  either  a  hit  or  fuel  exhaustion* 
Assumption  2:  At  the  time  of  acquisition,  the  target  is  opening  at 
maximum  speed.  Target  course  and  speed  are  constant  throughout  the 
attack  and  reattack. 

Assumption  3:  On  the  attack  and  after  reacquisition  on  each  reattack 
the  weapon  maintains  contact  and  closes  at  attack  (maximum)  speed* 
Assumption  k:  The  target  operates  in  depth  within  specified  minimum 
and  maximum  limits  such  that  the  maximum  depth  limit  is  not  greater 
than  the  weapon  maximum  operating  depth* 

Assumption  5:  The  depth  at  which  a  hit  may  take  place  is  a  random 
variable  denoted  the  hit  depth  which  has  an  assumed  probability 
density  function. 

Assumption  6:  The  average  weapon  depth  for  the  attack  is  the  arithmetic 
average  of  the  search  depth  and  the  average  hit  depth*  The  weapon  and 
target  depths  for  the  final  stage  of  the  attack  and  throughout  each 
reattack  are  equal  at  some  depth  within  the  limits  of  assumption  four* 
Assumption  7:  The  attack  and  each  reattack  is  from  astern  the  target* 
There  is  some  close-in  range  at  which  the  weapon  sensors  lose  contact 
with  the  target.  On  losing  contact  with  the  target  the  weapon  continues 
on  a  straight  track  at  attack  speed  for  a  preset  dead  time*  Upon  com- 
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pletion  of  the  dead  time  run,  the  weapon  instantaneously  slows  to 
search  speed  and  commences  a  left  turn  at  reattack  turn  rate  and 
constant  depth. 

Assumption8 :  The  weapon  reattack  turn  is  a  circle  with  no  advance  or 
transfer.  The  weapon  reacquires  the  target  upon  completion  of  one 
complete  reattack  turn  and  instantaneously  increases  to  attack  speed. 
Assumption  9:  At  the  time  the  close-in  weapon  sensors  lose  contact 
with  the  target,  there  exists  an  intersection  of  the  weapon  track 
extension  with  the  vertical  plane  through  the  submarine  target  center 
which  is  perpendicular  to  any  plane  containing  both  the  target  and  the 
weapon.  Let  the  target  projection  in  the  vertical  plane  be  approximated 
by  a  circle  denoted  the  target  circle  centered  on  the  target  longitudinal 
axis  with  a  diameter  equal  to  the  length  of  the  beam  dimension  of  the 
target.  Let  the  circle  center  be  the  origin  of  a  doubly  infinite 
axis  containing  the  intersection.  Assume  that  the  intersection  co¬ 
ordinate  on  the  axis  is  a  normally  distributed  random  variable  with 
mean  zero  and  standard  deviation  o  «> 

Assumption  10:  A  hit  occurs  if  the  weapon  track  at  the  closest  point 
of  approach  to  the  target  center  passes  within  a  distance  of  the  center 
not  greater  than  the  radius  of  the  target  circle  plus  a  small  proximity 
correction  factor. 

Define : 
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The  derivation  is  analogous  to  that  for  acquisition  probability 0 
Define: 


t  =  time  the  weapon  sensors  lose  contact  with  the  target  on  an 
attack  or  reattack  (zero  reference  time) 

RH  =  range  from  the  weapon  to  the  target  center  at  time  tQ  (yd) 
t  =  time  from  t  of  closest  point  of  approach  of  weapon  and 

L»iA  O 

target  (sec) 

MRH  =  distance  traveled  by  weapon  from  time  tQ  to  time  t^p  ( yd ) 

S  =  weapon  attack  speed  (yd/sec) 

WA 

S_,  =  target  maximum  speed  (yd/sec) 

BW  =  target  beam  width  (yd) 

VT  =  weapon  proximity  correction  factor  (yd) 

r  =  coordinate  of  intersection  of  the  weapon  track  extension 

with  the  vertical  plane  (A)  through  the  submarine  target  center 
perpendicular  to  any  plane  containing  both  target  and  weapon  at 
time  tQ  (yd) 

r*  =  coordinate  of  intersection  of  the  weapon  with  the  vertical  plane 
(A*  )  parallel  to  plane  A  containing  both  the  submarine  target 
center  and  weapon  at  time  t^p  (yd) 

F^(x)  =  cumulative  normal  probability  distribution  with  mean  zero  and 
standard  deviation  one  evaluated  at  the  point  x 
Consider  the  attack  situation  depicted  in  Figure  17« 

The  xy  plane  is  the  vertical  containing  the  weapon  and  target  at  time 
tQ*  Planes  A  and  Af  are  parallel  to  the  yz  plane  and  contain  the 
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Figure  17 %  Planes  A  and  A8  in  the  attack  situation 


bases  of  two  right  circular  cones  with  a  common  apex  at  the  origin 0  The 
line  OR  is  the  weapon  track  and  its  orientation  with  respect  to  the  x 
axis  determines  the  orientation  of  the  axes  LL  and  L8  L0  on  which 
coordinates  r  and  r1  are  measured.  The  common  intersection  of  the  cone 
of  base  radius  r*  9  the  plane  A1 9  and  the  xy  plane  is  a  line  of  length 

2 | r 1 1  .  By  similar  triangles 

i  MRH 
r  =  RH  r 

Since 


MRH  =  RH  S„A  /  (S„A  -  STM) 


then 


r'  =  r  SWA  /  ^SWA  "  STM^ 


Thus  r1  is  a  random  variable  with  mean  zero  and  standard  deviation 


“  0  SWA  /  ^SWA  “  STM^ 


Therefore 


h  = 


P  2  |  r*  |  <  BW  +  2  Vt] 

=  P  |r'  |  £  (BW/2)  +  VT] 


=  2  F 


0 


(BW/2)  +  VT 
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(BW/2)  +  VT 
0  SWA  /  ^SWA  "  STM^ 


Define: 

HD  =  hit  depth  (ft) 

^HD^X^  =  probability  density  function  of  the  random 
variable  hit  depth 
D^  =  maximum  target  depth  (ft) 

DMIN  =  minimum  target  depth  to  be  considered  (ft) 

Define  a  set  of  Q  equal  depth  subintervals  in  the  depth  interval 
DMIN*  °TmJ  suc^  fchat  the  depth  at  the  lower  boundary  of  the  q  th 

°MIN  +  q  DTM~DMINj‘ 

=  P  hit  depth  is  in  the  q  th  depth  subintervafj  j  q-1^  ° ° ° >  Q 
Assume  as  in  Figure  18: 


subinterval  is 


Define 


~  2  (x  -  Dw—r)  /  (D- 


fcHDw  -  -  -  -MIN 7  '  VUTM  "  Dimy 


Then; 


S3  -  l. 
Q2  ’ 


q 


i. 


Q 


Define : 

w'  =  fuel  used  in  the  n  th  attack  or  reattack  (lb) 
2n  ' 

=s  initial  acquisition  range  (yd) 
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fHD^x' 

A 


2 

dtm"dmin 


0 


Figure  18: 


Assumed  probability  density  function  of  the 
random  variable  hit  depth. 
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w^t  =  weapon  fuel  expenditure  rate  at  attack  speed  and 
at  the  average  depth  for  the  attack,  (ib/sec) 

HD  =  average  hit  depth  (ft) 

Then 

5B  =  2/3  Dra  ♦  1/3  dmin 

and  the  average  depth  for  the  attack  is 


WS 


+  HD 


-  1/2  Dwg  +  1/3  Dtm  +  1/6  D. 


TM 


MIN 


For  the  attack 


w 


W  = 
21 


2t 


WA 


TM 


Define: 


w 


Stq 


w 


Atq 


'DT 


'CT 


=  weapon  fuel  expenditure  rate  at  search  speed  and  at 
the  midpoint  of  the  q  th  depth  interval .  (lb/sec) 

=  weapon  fuel  expenditure  rate  at  attack  speed  and 
at  the  midpoint  of  the  q  th  depth  interval,  (lb/sec) 

=  dead  time  from  t  at  which  the  weapon  commences  the 
reattack  turn  (sec) 

=  time  from  t  at  which  the  weapon  completes  the  reattack 
turn  and  reacquires  the  target  (sec) 

=  reattack  rate  of  turn  (deg/sec) 


Then; 


t 


CT 


'DT 


+  360° 


t 


CPA 


RH 

SWA  '  STM 
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The  weapon  will  spend  the  following  time  at  attack  speed  during  each 
reattack:  See  Figure  19. 

^TCT  "  fcCPA^  STM~  ^  tDT°tCPA  \  SWA 
A  =  DT  -  CPA  +  SWA  - 

The  weapon  will  spend  the  following  time  at  search  speed  during 
each  reattack: 

T  -  360° 

S  ^ 

Therefore : 


Q 

W2n  *  21  [TA  WAtq  +  TS  “stq  ]  PEq 
q=l 


n  =  2,  •••,  N 
Define : 

S  =  minimum  preset  value  of  attack  speed  (yd/sec) 
o 

w^  =  available  fuel  weight  for  hit,  event  2,  (lb) 

w”  =  available  weight  for  modifications  to  the  propulsion 

system  to  provide  a  maximum  speed  capability  in  excess 

of  S  (lb) 
o 

N  is  a  function  of  w^  according  to 


N 


N+l 


w_ 


2n 


—  << 


w, 


2n 


n=0 


n=0 


However,  recall  that  w^n  is  a  function  of  attack  speed  which  is  it= 
self  a  function  of  w^' . 
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Figure  19:  Geographic  plot  of  weapon  and  target  tracks 

separated  to  show  time  differences •  The 
distancesbelow  the  target  track  are  measured 
from  the  joint  position  at  time  t  . 

v  iA 
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Thus 


N(v£,  w£'  ) 

Pr(N)  =  1-  [l-hn^') 

n=l  L. 


Recall  w2  =  +  w^' 

I  _l 

W,  W  =  lower  and  upper  limits  respectively  on  the  range  of  fuel 
weight  (lb) 

m  Jl 

W  j  W  =  lower  and  upper  limits  respectively  on  the  range  of 
weight  for  modifications  to  the  propulsion  system  to 
provide  a  maximum  speed  capability  in  excess  of  Sq  (lb) 
For  each  value  in  the  range  for  w^ 
min  -£w'  ,  W*  'j-  <w2<~min  W  +  W  ,  v/}- 

the  function  P  is  evaluated  for  each  possible  combination 
H 

of  the  variables  w^  and  w^f  as  they  are  allowed  to  vary  such  that 

they  remain  within  their  respective  ranges 

1  1 

W  ^ w.  ^  W 
-  —  2 

it  ii  _i  i 

W  <w2^W 

and  their  sum  is  equal  to  w2*  Let  K  be  the  number  of  values  of 
w2«  This  procedure  provides  sets  of  values 

w*  M  }  w’  "  .  P„  M  y  j  =  l,  *  *  * ,  J, 

2j*  2ik*  Hjk'  J  7  7  k 

for  the  k  th  value  of  w^*  (k  =  1,  • 0  •  9  K) 

The  entire  set  of  sets  of  values 


W2k*  W2jk'  W2jk>  PHjk;  j  =  1# 
k  =  1,  •**,  K 


O  Q  O 

> 


J 


k 
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is  the  input  to  the  dynamic  programming  allocation  model 0 

Input  Requirements 

1.  Target 

BW  =  target  beam  width  (yd) 

DTM  =  target  maximum  depth  (ft) 

=  target  maximum  speed  (yd/sec) 

2 .  Weapon 

S  =  weapon  search  speed  (yd/sec) 
w  j 

VT  =  weapon  proximity  correction  factor  (yd) 

00^  =  weapon  reattack  turn  rate  (deg/sec) 

i  _i 

W,  W  =  lower  and  upper  limits  respectively  on  the  allowable 
range  of  fuel  weight  (lb) 

Sq  =  minimum  allowable  weapon  attack  speed  (yd/sec) 

ii  _i  i 

W  ,  W  =  lower  and  upper  limits  respectively  on  the  range 
of  weight  for  modifications  to  the  propulsion 
system  to  provide  a  maximum  speed  capability  in 
excess  of  S  (lb) 

G^  =  Relationship  between  increased  weapon  speed  above  Sq  and 
increased  weapon  weight 

=  Relationship  between  fuel  expenditure  rate  and  speed  and 
depth.  Recall  that  this  relationship  may  vary  as  a 
function  of  attack  speed  only  in  the  range  for  speeds 
above  search  speed.  This  relationship  is  required  to 
specify  and  w2j.  (lb/sec) 
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3.  Miscellaneous 


RH  =  range  from  the  weapon  to  the  target  center  at  the  instant 


(tQ)  during  close-in  attack  when  the  weapon  sensors  lose 


contact  with  the  target  (yd) 

0  =  standard  deviation  of  the  normally  distributed  random 

variable:  coordinate  of  intersection  of  the  weapon  track 
extension  with  the  vertical  plane  (A)  through  the 
submarine  target  center  perpendicular  to  any  plane  con¬ 
taining  both  target  and  weapon  at  time  t  •  (yd) 

DllTKT  =  minimum  target  depth  to  be  considered  (ft) 

MIN 

=  initial  acquisition  range  (yd) 

Q  =  number  of  depth  subintervals  in  the  depth  interval 


Intelligence  reports  and  estimates  or  CNO  requirements  may 
provide  the  target  inputs.  The  weapon  inputs  are  available  from 
weapon  designers  at  least  in  the  form  of  educated  estimates.  The 
inputs  RG  and  O  are  functions  of  the  medium,  the  weapon  sonar  and 
guidance  systems,  the  target  physical  dimensions,  and  the  orientation 
relationships  between  the  target  and  weapon.  In  the  absence  of  a 
specific  theory  these  must  be  handled  by  simulation  or  sensitivity 
analysis.  The  input  is  a  function  of  the  capability  of  the 

weapon  delivery  system  and  the  tactical  situation  as  seen  by  the 
target  commander.  In  the  absence  of  specific  considerations  cri¬ 
terion  one  would  determine  to  be  zero.  The  input  RA  should  be 

MIN  A 
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selected  close  to  the  maximum  weapon  capability  in  accordance  with 

criterion  one*  The  function  P  may  be  very  sensitive  to  changes  in 

R  and  an  examination  of  results  over  a  range  of  values  for  R  may  be 
A  A 

appropriate*  The  input  Q  should  be  large  enough  to  point  up  the 
effect  of  an  irregularly  shaped  fuel  expenditure  rate  curve  but  small 
enough  not  to  excessively  prolong  the  computer  program* 
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CHAPTER  V 


PROBABILITY  OF  KILL  DAMAGE 

Consider  a  tactical  situation  in  which  a  weapon  explodes  within 
a  specified  proximity  of  some  portion  of  a  submarine  target  and  make 
the  following  assumptions* 

Assumption  1:  The  hit  depth  is  a  random  variable  which  has  an 
assumed  probability  density  function* 

Assumption  2:  The  point  on  the  target  at  which  the  hit  occurs  de¬ 
noted  the  hit  point  is  a  random  variable  which  has  a  uniform  distri¬ 
bution  over  the  entire  surface  area  of  the  target*  The  weapon 
approaches  the  target  from  astern,  loses  contact  at  some  close-in 
range  and  continues  on  a  straight  track*  Since  the  target  may  be 
in  a  turn,  climb  or  dive,  all  points  of  the  target  are  available 
as  hit  points  * 

Assumption  3:  The  target  operates  in  depth  within  specified  minimum 
and  maximum  limits  such  that  the  maximum  depth  limit  is  not  greater 
than  the  weapon  maximum  operating  depth*  Divide  the  target  surface 
area  into  P  area  elements,  each  of  which  has  a  constant  hull 
vulnerability  throughout.  Divide  the  interval  of  possible  hit 

depths  into  Q  equal  depth  subintervals  such  that  the  depth  at  the 

i 

lower  boundary  of  the  q  th  subinterval  is 


dmin  +  q 


dtm"dmin 


82 


where 


DTM  =  maximum  target  depth  (ft) 


MIN 


=  minimum  target  depth  to  be  considered  (ft) 


Define:  (w^)  =  P  £kill  damage  given  an  explosive  weight  wjj 

f„  (x)  =  probability  density  function  of  the  random 
HD 

variable  hit  depth 

P£^  =  P  £hit  depth  is  in  the  q  th  depth  subinterval J  | 

q  =  i>  Q 

P^  =  P  j^hit  point  is  in  the  p  th  surface  area  elementj; 


P  =  1> 


P 


=  hull  vulnerability  for  kill  damage  of  the  p  th 
target  surface  area  element  at  the  midpoint  depth 
of  the  q  th  depth  subinterval  given  an  explosive 
weight  w^.  This  number  is  one  if  kill  damage  occurs 
under  the  specified  conditions  and  is  zero  if  kill 
damage  does  not  occur • 
w^  =  explosive  weight  (lb) 

It!  _ V  «  • 

W  ,  W  =  lower  and  upper  limits  respectively  on  the  range 
of  explosive  weight  (lb) 

Assuming  as  in  Chapter  IV 


£hd(x)  - 


2  (*-Dmtl) 

(dtm*  dmin>! 
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Then 


P  =  —  *1  — - •  a  =  1  •••  o 

%  q2  '  q  h  ’  * 

p  _  area  of  p  th  surface  area  element;  p  s  1,  P 

Ep  “  total  target  surface  area 


PE  =  PE  ^w3  ^  = 


HV  (w.)  P_  P_ 
pq  3  Ep  Eq 


p=i  q=i 


I  M  M 


The  table  of  P  versus  w  as  w  varies  from  W  to  W 
b  5  5  — 

is  the  input  to  the  dynamic  programming  model* 


Input  Requirements 

1.  Target 

PEp  =  P  j^hit  point  is  in  the  p  th  surface  area  element] 
p  =  1,  •  •  •,  P 


=  maximum  target  depth  (ft) 

TM 

2.  Weapon 

iii  _i  i  • 

W  ,  W  =  lower  and  upper  limits  respectively  on  the  range 
of  explosive  weight  (lb) 

3o  Target  and  Weapon 

HV  (w^)  =  hull  vulnerability  for  kill  damage  of  the  p  th 

pq  3' 

target  area  element  at  the  midpoint  of  the  q  th 
depth  subinterval  given  an  explosive  weight  (w^) 

4*  Miscellaneous 

D  =  minimum  target  depth  to  be  considered  (ft) 

MIN 

Intelligence  or  CNO  requirements  may  provide  D_w.  P  and  some 

TM  Ep 

of  the  information  for  HV  (w^)  .  Weapon  designers  may  estimate 

it*  _i i i  pq  j 

W  and  W  and  weapons  engineers  may  provide  extrapolated 
estimates  of  HV  (w_). 

pq  3 
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APPENDIX  A 


FLOW  CHART  NOTATION 


The  following  notation  is  used  in  the  flow  charts  of  the  remaining 
appendices . 


Function  Box 


A  — 


The  function  box  contains  a  number  or  expression  to  be  evaluated 
A  which  is  to  be  placed  in  memory  cell  Bo 

Choice  Box 


The  choice  box  contains  a  comparison  of  the  number  or  express 
sion  A  to  the  number  or  expression  B,  Exits  are  chosen  as  follows s 

Exit  is  chosen  if 
A>  B 
A  >  B 
A  =  B 
A  <B 
A  <  B 


Exit  Label 


< 

< 
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In  Box 


INPUT 


The  in  box  contains  an  input  to  the  program  *  In  appendices 
C  1  through  C  8  the  output  from  one  appendix  is  carried  over  intact 
into  the  next  appendix.  Inputs  to  be  used  for  the  first  time  in  the 
sequence  of  flow  charts  are  placed  in  a  separate  in  box* 

Out  Box 

OUTPUT 


The  out  box,  one  for  each  appendix,  contains  the  output  of  the 
appendix  program. 


Index  Box 


{ 


1  +  i 


} 


The  index  box  contains  a  manipulation  on  an  integer  in  an 
index  register  and  is  the  tool  which  provides  for  loops  in  the 
program . 
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Subroutine  Box 


The  subroutine  box  indicates  a  return  jump  to  the  subroutine  as 
numbered  by  k.  The  flow  chart  of  each  subroutine  is  shown  after 
the  flow  chart  of  the  first  routine  using  it , 

Return  Box 


A  return  box  indicates  the  end  of  a  subroutine  implying  a  jump 
back  to  the  main  routine. 

Fixed  Connector 


A  fixed  connector  is  used  to  connect  remote  parts  of  the  flow 
chart  with  one  another. 


Assertion  Box 


Note  j 


- 1 

An  assertion  box  refers  the  reader  to  a  note  found  at  the  end  of 


the  flow  chart. 
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Instruction  flow  is  generally  from  left  to  right  and  from  top  to 
bottom.  All  boxes  are  entered  from  the  left  or  from  the  top  and  all 
exits  from  boxes  are  made  down  or  to  the  right.  The  arrow  is  used  to 
further  indicate  the  sequence  of  computer  instructions.  Arrows  are 
used  at  all  line  intersections  to  indicate  entry  into  the  program 
between  two  boxes. 
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APPENDIX  B 

THE  ALLOCATION  MDDEL 
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Note  1:  Provision  is  made  for  a  set  of  payload  weight  limits  (W  )  in 

n 

case  it  is  desired  to  examine  a  range  of  values.  The  input 

W,  is  to  be  the  greatest  of  the  W  . 

N  n 

Note  2:  The  subroutine  RJl  computes  the  h^(w)  and  the  w„(w)  for 
i=l,  3.  h^  is  handled  separately. 
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Note  3 


Note  Ij- 

Note  5 
Note  6 

Note  7 
Note  8 

Note  9 


This  loop  provides  that,  for  the  W  considered,  the  function  tr 

is  maximized  over  the  w^ .  For  h^  the  maximum  meaningful  value 

in  both  loops  is  w^=D  so  the  subroutine  is  stopped  by  the 

comparison  W  :  D.  However  for  h^,  the  maximum  meaningful 

value  in  the  two  loops  is  different,  the  one  in  this  loop  being 

w-._  and  the  one  in  the  other  loop  being  W...„  s  D. 

3  b  MAX 

This  loop  stores  the  h^(w)  and  w^(w)  obtained  as  W  is  allowed 
to  vary  over  its  range  of  meaningful  values, 

This  loop  maximizes  the  w^^  over  j  for  a  given  k„ 

The  maximum  meaningful  value  of  w^  for  a  given  W  is 


min 


The  h^(w),  w^Cw),  and  w^’Cw)  are  stored  for  each  value  of  W. 

This  loop  works  backwards  through  the  stages  to  retrieve 

the  optimum  weight  values  for  each  input  payload  maximum 0 

The  are  the  optimum  volume  sums  to  be  compared  with  the 

input  payload  maximum  allowed  volume  V0  If  the  solution  is 

not  optimal  with  respect  to  the  restraints^  a  new  X  is  chosen 

and  the  program  repeated.  If  the  final  solution  is  obtained^ 

then  fuel  weight  is  +  1^,  weight  for  added  attack  speed  is 

H  and  explosive  weight  is  P  « 
n  n 
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APPENDIX  Cl 

AVAILABLE  ECHO  EXCESS  (p^) 
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Note  1;  This  flow  chart  neglects  the  effects  of  ambient  sea  noise  and 
of  volume  reverberation.  Both  are  considered  to  have  negli¬ 
gible  effect  on  the  interference  background  masking  level. 

If  it  is  desired  to  consider  the  volume  reverberation,  it  is 
M(Rv)k  =  WS  -  21^  +  10  log  mv  “(Jv)k  +  10  log 

+  20  log  (100k)  +  RD;  k  1,  •«»,  K 

where  M(R  )  =  volume  reverberation  masking  level  at  range  100k  yards  (db) 

'  v  k 
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WS  =  on~axis  source  level  (db  at  1  yard  from  transducer) 

=  one  way  transmission  loss  at  range  100k  yards  (db) 

10  log  m^  =  volume  scattering  coefficient  (db) 

(j  )^  =  volume  reverberation  index  of  the  transducer  (db) 
c  =  velocity  of  sound  in  the  medium  (yd/sec) 

AT  =  ping  length  (sec) 

RD  =  recognition  differential  of  transducer  (db) 


The  index  (jy)^  may  be  approximated  by 


=  200  ko) 


(Jv)k  =  "10  los  b3  \  (  P  -s }d  _n_ 


k  =  l. 


K 


where  m  =  search  rate  of  turn  (deg/sec) 

j  =  off-axis  angle  in  the  vertical  plane  (deg) 

P  =  off-axis  angle  in  the  azimuth  plane  at  time  of 
pulse  transmission  (deg) 

b^,  ^2*  s  normalized  pattern  functions  of 

the  transmitting  (l,  3)  and  receiving 
(2<>  4)  patterns  in  the  vertical  (1,2) 
and  azimuth  (3$  4)  planes 

_Q_  =  solid  angle  subtended  by  the  sound  energy  beam0 

This  in  effect  assumes  a  medium  of  infinite  extent 
in  the  direction  of  the  sound  beam0 
Note  2s  M(SN)  =  SN  +  10  log  A  f-DR  +  RD 

where  M(SN)  =  self-noise  masking  level  (db) 

SN  =  equivalent  isotropic  self-noise  level  at  search  speed 
and  depth  (db) 
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Note  3t 
where 


Note  kt 
where 

Note  5: 

Note  6: 
where 

Note  7: 
where 

Note  8: 


£  =  receiving  bandwidth  (cycles/sec) 
DR  =  receiving  directivity  index  (db) 
2H  «  WS  +  TS-M(SN) 


H  =  maximum  possible  on^axis  range  for  acquisition  (yd). 

This  neglects  noise  levels  except  the  self=noise  level 
and  neglects  the  weapon  enabling  level. 

TS  =  target  strength  (db) 

=  20  log  100k  +  100  ak 

a  -  attenuation  coefficient  (db/yard)  and  it  is  assumed  that 
the  transmission  anomaly  is  zero. 

This  loop  finds  the  upper  limit  (lOOK)  on  the  acquisition 
ranges  to  be  considered. 

M(TN )k  =  TN  +  10  log  A  f-Hk  +  RDj  k-1,  •••,  K 
M(TN)k  =  target  noise  masking  level  (db) 

TN  =  radiated  noise  level  of  the  target  (db  at  1  yard  from 
transducer) 


D. 


<#> .  =  arc  sin 
k 


WS 


300k 


=  the  angle  above  the  horizontal  at  which  an  arc  of  radius 

100k  yards  and  of  origin  at  search  depth  will  inter= 

sect  the  surface  of  the  water 

D  „  =  weapon  search  depth  (ft) 

W  o 

(j  )  =  surface  reverberation  index  of  the  transducer  (db) 

\i  “angle  of  attack  (pitch)  of  the  weapon  (deg) 

This  is  an  approximation  described  in  Chapter  k-  of  refer- 
ence  k . 
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Note  9:  M(rs)^  =  surface  reverberation  masking  level  at  a  range  of 

100k  yards  (db) 

10  log  m^  =  surface  reverberation  scattering  coefficient  (db) 
Note  10:  This  subroutine  adds  the  effects  of  two  noise  levels  in 

decibels.  A  rounding  off  process  is  employed  based  on  Figure  3 
of  reference  3*  A  difference  in  levels  of  more  than  10  db  pro¬ 
vides  a  sum  effect  equal  to  the  larger  of  the  individual 
levels.  Two  equal  levels  have  a  sum  effect  equal  to  3  db 
higher  than  either  individual  level. 

Note  11;  WEL^  =  weapon  enabling  level  (db).  Values  out  to  a  range 

of  100K1  yards  are  inputs  where  is  estimated  such 
that  it  will  be  larger  than  the  K  found  in  the  program. 
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APPENDIX  C2 

SCANNED  VOLUME  RADIALSp  GENERAL  CASE 


INPUT 


INPUT 


PRRfl ,  « 9  « '  p  D, 


TM 


APPENDIX  Cl  OUTPUT  j - aytRR^-^b 


200  a^/c  { 1  >  1  y~ =  Kq  -c-  (i  =>  1)  0^ — 


P 


Note  1 


-(  1  ♦  i-^i> 


T> 


i  -  1^1 


>  ~»J> 


f 


Note  2 


RJ3) - <l->k 


p  j2k  0  b  + 


/  1  »  k  ->¥V 


r1 — -  (iib^n  ♦  ♦  a 
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RJ4) 


1  — in  -  '  1  ^  k  )  —  ^  0  ja  x - 1 


i 


pjmk  :  pJm(k+l) 


j'Note  9 1 


(~1  ±  k  ->k  ) — (  1  »  n  — 5>n  / — 


P 


Note  10 


Pjak-^PJ(k-n)hL-  |100k  •  c/2PRRa  p^  c/2PRRs-^Rj(k^} 


100k  — ->■  R 


'J  (k-n) 


J(k=*n)  •  -sin(j+n)ees 


(dtm  -  Dws}/3 


R, 


sin  ( j  +ji)cos^k^D^  J  (k=m) 


-  V* 

p  • 

> 

“u/3 

^JUc-n)  sin(J+|i)cosp^ktta^ 

sin(J^)cospj  ~^Rj  (k»n) 

[Note  12 


pj(k-n)  »  p’ 

< 

(TTH 

:^k) 

k  -  n  — ^K„ 


(RTN 


Pjl(k-n)  *  pJ2k 


<  k  -  1  -5>k^ - ( n  2  — >a) - (J) 


-  1  — $»n) - (2 -g» m  )- 


Note  13 


>  -  l^»k) - <2_^m) - (n_-2_- ^n> 


<D 
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which  meet  the  requirement  that  the  echo  will  be  received 
in  the  interval  I  —  <*9  j  • 


Note  2:  The  overlap  in  plane  zero  is  computed  first  to  provide  the 
angle  p*  .  Then  p*  <=  b  and  p®  define  the  limiting  angles 

of  the  P  jmk" 

Note  3:  AJki  ■  Ajk(  P  )  for  the  i  th  P  . 
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Note  4; 

Note  5; 

Note  6s 

Note  7t 

Note  8; 

Note 

Note  10; 
Note  11; 


After  the  computation  of  which  indicates  that  the  maximum 

acquisition  range  in  plane  j  and  all  of  the  p  in  plane  j 

jmR 

have  been  computed^  the  subroutine  ends0 

This  indicates  that  there  is  no  acquisition  probability  in  the 
plane  j  and  changes  the  range  of  planes  to  be  considered  from 
-J*  °  •  •  9  J  to  j\  •  °  °  9  J5  5  . 

In  this  case  the  p  is  limited  by  the  P  c  and  would  have 
been  larger  if  the  overlap  between  pings  had  not  been  removed 
from  the  considered  scanned  area 0 

This  loop  starts  with  the  leftmost  angle  (  p 0  ®(l®l)  ^  ) 

to  be  considered  and  computes  each  angle  to  be  considered 

indexed  by  i  and  incremented  by  0^  9  to  a  maximum  angle  pTsP°  © 

This  point  is  the  commencement  of  the  loop  which  will  compute 

the  p  and  R  starting  with  j  *  ~J0 
jK  JK 

There  may  be  several  equal  P because  the  scan  azimuth  limits 
were  reduced  by  the  overlap 0  The  first  of  the  equal  p  ^  to  be 
considered  will  then  be  the  one  of  maximum  k0 

q 

The  ping  interval  range  is  yards/ping. 

The  range  at  which  the  radial  at  angle  p  ^  intersects  the 
surface  of  the  water  for  j  >•  -  jiis  computed  below. 

See  Figure  20.  The  distance  0C  is  desired.  The  search  depth 
is  (DWS/3)  yards. 


length  (CB)  =  Dyg  /  3  sin(j^p) 
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y 


Figure  20.  Computation  of  radial  length  for  angle  p 

The  plane  j  for  which  the  scan  outline  is  shown 

is  perpendicular  to  the  vertical  xy  plane.  Angles 

CBO  and  DBO  are  right  angles .  Length  Sis  tefee 

search  depth  of  Dyjg/3  yards. 

I 
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Note 


Note 


length  (OC)  =  length  (0B)/eospjk  =  Dyg/3  sin(jj-frji)  eosp 


For  j  -  the  development  is  the  same  except  that  (j  +  p*  ) 

is  negative  and  the  considered  depth  is  the  difference  between 

the  maximum  target  depth  and  the  weapon  search  depth 0 

TM 

12 s  If  to  the  left  of  the  radial  of  maximum  range  K.  then  m  -  1 

J 

and  the  following  angle  will  have  a  larger  kc 
If  to  the  right  of  the  radial  of  maximum  range  K.  then  m  ^  2 
and  the  following  angle  will  have  a  lesser  k0  There  will  be 
one  or  two  radials  of  maximum  range 0  That  is 


PjlKj  —  ^J2Kj 


13:  If  PjlK  =  Pj2R  »  the  program 


program  includes  it, 


8  o  Otherwise  the 

j2kj 
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APPENDIX  C3 


SCANNED  VOLUMES^  GENERAL  CASE 


Note  Is  SVq^  i-s  not  required  in  the  program  but  it  may  be  desired 
to  print  it  here  as  information. 


105 


APPENDIX  C4 


SPECIAL  VOLUME  ELEMENTS l  GENERAL  CASE 
INPUT  INPUT 


£H 


STS 9  SWD‘ 


^DS ! 


W» 


Ot9  wlt 


APPENDIX  C3  OUTPUT 


} 


^11360/^)11 


3 


+  1  ^  M 

*  -L  z**-  n 

PWSWOt^3  8ia  ^DS  SWD~>MOSD 


□ 


P'“W_n)/w. 


OSD  ""It 


PRRg/M]|  +  1^-U 


-TES 


]D- 


t£ 


»  0 


RJ6 


D« 


OUTPUT 

^du9  ^du9 
LBDV, 


du9 - dir  duv  9 

DVdu9  ^du9  RBDVdu9  d  *  19  ° 

vdu;  d  =  d+i9  ooo9  d»;  u  -  i9 
BM7duJ  d  =  D+l9  "\  D'5«=l, 

<9 

rr  o  j  —  i  o  o  ©  t\  o  ©  „ 


1»  00%  D®  u 


ASVr 


;  u  ^  29  oo 


L 


BB^0du9  d 

=  1, 

°°°9  D9 5 

u  : 

s  29 

°°°9  u 

qSV0uJ  q  = 

=  a9 

B9°  u  =  29 

0  0 

*9  U 

hAd9 

BAId 

9  tan^AId 9  d  = 

1. 

o  o  © 

9 

D* 

BAOd 

,  wAM; 

d  -  19  0 

0  0 

9 

D 

hBd9 

BBId 

9  9  d 

1, 

0  o  o 

9 

D« 

BBOd 

*  WBOd5 

d  -  19  0 

0  0 

9 

D 

c=> 

PJk9 

Rjk9 

j  = 

J«,  °°°9 

J»  « 

0  k  s 
9  K 

“ij  o  o  © 

x9  9 

^DS 9 

SWD9 

STS 9 

WOSD9  Wlt9 

m9  u 

P*» 

^9 

Dws9 

B 

j 
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L-^RTN) 
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Note  10 
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Note  13 


RJ13 — L- 


A2tan(6/2)  hq(J-^B 


(RTN 


RJ14H- 


I'Note  14 


(1/3)  tan (6/2)  C(D3  -  E3)^=>F 


RTN 


Note  11 


(RJ17) — 


(TES  =•  L(d=2))  =  BSd) 

(d-1)  “ 


HL 


HL, 


—  d-^D 


{YTT^rJ—  (HLd  -  BSd)/(Hl(d_1)  -  BSd)  ^»CT 


~ '  HL(d-l)  *  HLd  ^‘hqd 


(RTN 


depth. 

Note  2:  The  following  comparison  is  used  to  determine  the  plane 
j  to  be  used  as  a  reference  plane  for  the  computations. 


Ill 


Note  3t 


Note  k 


Note  5 


Note  6 


Note  7 


The  following  four  function  boxes  set  up  reference  constants 
for  the  subsequent  computations  <> 

:  The  index  d  is  used  to  identify  the  levels 0  The  index  (d«r) 
identifies  the  current  plane  j  above  the  plane  TJ«  The 
index  (d-s)  indentifies  the  current  plane  j  below  the  plane 
TJ-  ATJo 


i  This  comparison  determines  if  the  level  d  is  to  be  determined 
by  a  plane  j  >  TJ  or  by  a  plane  j  <(TJ=  AtL)0 

It  is  only  required  for  depths  at  or  above  the  critical  depths 
that  is  for  d  ^  D . 

:  This  subroutine  computes  the  distance  L^0  See  Figure  21 • 

As  shown  H^tj  +  ^  sin  (TJ  +  2  +  p )  > 

H(tj  g)  s^n  +  p )  and  ATJ  -  0.  The  distance  to  be 

computed  is 

L,  =  length  AB  =  length  FG  =  length  FE~  length  GE 
d 

Ld  =  length  FE  -  ift  (U,18th  HE) 

;  If  the  level  length  is  longer  than  the  distance  TES 

traveled  by  the  submarine  during  a  complete  weapon  search 

turn,  the  critical  depth  is  at  a  level  above  L,  but  below 

d 

L/  .  , ».  Then  the  level  defined  by  TES  must  be  found  and 
\d“l ; 

designated  and  the  original  level  dropped  from  con= 

sideration  until  later.  If  L  is  shorter  than  TES  then  it  is 

above  the  critical  depth  and  its  parameters  must  be  computed 

to  find  qSV  ,  ,  the  volume  searched  in  the  d  th  level  on  the  u  th 
odu' 

search  turn,  (q  =  A,B). 
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TJ 


Figure  21:  A  sample  volume  portion  for  computation  of  L , .  The  points 

d 

(TJ^  d)  are  the  radial  endpoints  of  the  corresponding 


planes.  The  level  is  under  consideration  and 


L3>  1s 

to 

be  computed. 

The 

bc3- 

SC3 

=  length 

FE 

BS3 

ss3 

=  length 

BG 

ss3 

ss3 

=  length 

HI 

SC3  " 

sc3 

-  length 

HE 

113 


Note  8:  See  Figure  22.  Subroutines  RJ8,  9/  and  10  compute  the  distances 

and  angles  on  the  right  portion  of  the  figure  which  have  not 

already  been  computed «  Then  subroutine  RJ11  computes  all  of 

the  special  volumes  in  level  After  RJ11  the  level  index 

as  well  as  the  index  r  or  s,  is  increased  by  L  If  index  s 

is  changed  it  means  that  the  plane  which  determined  the  level 

was  above  plane  TJ.  Thus  d  and  s  must  be  increased  together 

so  that  the  difference  is  constant.  It  is  their  difference 

that  is  used  for  the  comparison  described  by  note  5° 

Note  9 l  This  subroutine  computes  the  scanned  volume  in  level 

(where  d^D)  on  each  search  turn  after  the  first  .  The  first 

comparison  is  used  to  find  out  if  the  u  th  search  turn  is  the 

last  effective  search  turn  in  level  L^0  If  it  is  not*  then 

subroutines  RJ12,  15/  and  16  are  used  to  compute  the  volumes 0 

The  cross-section  in  this  case  looks  like  that  of  the  right 

portion  of  Figure  22.  If  the  u  th  search  turn  is  the  last  of 

the  effective  search  turns  in  level  L,,  then  the  normal  cross® 

a 

section  is  cut  at  some  point  as  in  the  left  portion  of  Figure  22o 
In  this  case,  the  point  where  the  cross®section  is  cut  is 
found  by  one  of  the  subsequent  three  comparisons  and  the  cal¬ 
culation  of  scanned  volume  is  modified  accordingly 0 
Note  10:  This  corresponds  to  a  vertical  through  the  transducer  which 

goes  through  a  point  in  level  such  as  point  S  in  Figure  22 „ 
Note  11:  This  is  used  to  compute  the  height  of  the  critical  depth 

above  the  search  depth  (length  OR  in  Figure  22)  and  a  correc® 
tion  factor  to  be  used  to  find  the  constants  for  level 


114 


Figure  22:  Vertical  cross-section  through  levels  Lg  and  L ^ . 

The  right  portion  depicts  the  volume  scanned  on  the  first 
search  turn.  The  vertical  OAB  extends  down  to  the  trans¬ 
ducer  . 

Volumes  corresponding  to  depicted  areas  are 
Area  Volume 


HID 

DEIJ 

EJK 

CFG  =  EJK 
ABCF  ^  DEIJ 
HDEK 
ABCG 

The  distances  used  to 
length  NQ  =  TES 
length  BD  = 
length  BE  =  Bqo(d_1} 
length  AH  =  Bq^d 
length  AK  =  Bqod 
length  HK  =  Ld 


LqDV 

MqDV 


dl 


dl 


RqDV 


dl 


R,DVdy  i  RqDVdl 
MqDVdu  ±  MqDVdl 


qsv 

qSV 


odl 


odu 

compute  volumes  in  level  L,  are 

d 


length  HI  =  b 


Id 


length  JK  --  length  FG  -  b 


od 


length  AB  s  hqd 
length  AR  ^  HL^ 
length  BR  =  HL^ 
length  DE  -  ^ 
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Note  12:  This  is  the  volume  computation  for  FqDV^.  This  computation 
begins  by  turning  the  volume  upside  down  so  that  level  L,  is 

U 

above  level  Then  the  computation  is  the  same  as  that 

for  (corresponding  to  area  ABDH  in  Figure  22)  required 

later . 

Note  13:  See  Figure  23. 

Note  lh:  See  Figure  2b. 

length  A  I  =  D  tan# 
area  IJK  =  D2  tan  (6/2) 

volume  AJKI  =  (l/3)  (area  IJK)  (length  Al) 

=  (l/3)  tan  (6/2)  tan# 


Note  15: 


Note  16: 


Note  17: 


Note  18: 


The  volume  computations  (qMV^)  for  levels  below  the  critical 
depth  are  begun  at  this  point. 

The  volume  computations  for  the  last  level  use  B  for  the 

previous  level  and  the  angle  j  +  d  ■=  r  +  p  because  there  is 

no  plane  j  on  which  to  base  the  calculations. 

The  calculations  below  search  depth  are  identical  to  those 

above  after  the  H.  below  are  substituted  for  the  H.  above. 

J  J 

The  final  section  computes  the  scanned  volume  on  the  u  th 
search  turn  separately  for  above  and  below  search  depth. 
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Figure  23:  Volume  computation  of  a  triangular 


EB 


=  A.  Length  BC  =  A  tan  (6/ 2).  Length  ED  ■  h^ 
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w 

•H 
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The  altitude  lengths  are  AI  and  AH 


APPENDIX  C5 


FIRST  TURN  SCANNED  VOLUME:  i  TH  PING 


Note  1:  p  is  now  the  angle  between  Rq  and  the  radial 
(Rjk, 
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APPEND IX  C6 


SCANNED  VOLUMES  SUBSEQUENT  TO  FIRST  TURN:  i  TH  PING 
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(?) - |TId  -~fa5jTEsTCT 


Note  7 


> 

( Note  8 

1 - - 

1 

0 


(CH)2tan(&/2)  [hqd  -  (BqI(d-1j  -  (u-l)TES  -  CH)tan0qIdj 
-  (l/3)tan(&/2)W>qia  [(CH)3  -  (BqJfl  -  (u-l)TES)^] 


+  ASV— ASV 


(?) — (  1  d  d  ) - •  •  (?) 


Note  9 


|B_w,  r  (u-2)TES  l  CH 


Note  10 


n> 

1 

L  BqIa  -  (u-2)TC?Tch] 

i<J 

(Note  11 

1  ^ 

*ST  1SToa„^S7l - O 


(l/3)tan(&/2)tan£qId  [(CH)3  -  (B  H  -  (u-2)TES)3] 

-  qM?du  +  ASV  -5*  ASV 


'Note  11 

l> 

□ 


(CH)2tan(6/2)h  d  -  qM7du  *  ASV  ASV 


|BqId  -  (u-l)TES  :  CH 
”> 


,< 

'Note  12 

I' 

1 

0 


(CH)2tan(b/2)[hqd  -  (Bql(a_1)  -  (u-l)TES  -  CH)tar>0qld] 

-  (l/3)t«n(6/2)tan^qla  [(OH)3  -  (BgH  -  (u-l)TES)3]  *  tSV-3-iSV 

3=5) 


©-C 


1  +  d  -5»-  d 


L|  ASV  ^qSVlu  ) - (  1  +  u  — >•  u 


.  f 


Note  14 


j'  Note  13 


0 


q  :  A>=_L 


"© 


<B^q 


D  -^.d. 


D«  -5»-d, 


0 
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Note  1 
Note  2 

Note  3 


Note  4 
Note  5 


Note  6 

Note  7 

Note  8 


P  is  now  the  angle  between  Rq  and  the  bisector  of  the  ith  ping* 

:  In  this  case  the  target  volume  cylinder  is  large  enough  such 

that  the  scanned  volumes  after  the  first  search  turn  are  not 
effected * 

i  In  this  case  the  target  volume  cylinder  is  small  enough  such 

that  some  of  the  scanned  volumes  after  the  first  search  turn 

must  be  modified *  This  is  accomplished  above  and  below  the 
search  depth  and  for  each  search  turn  by  summing  the  modified  or 
unmodified  volumes  over  the  levels  to  obtain  qSV^*  First 
each  search  turn  above  search  depth  is  considered*  then  each 
search  turn  below  * 

In  this  case  the  original  qSVQ^  is  not  changed <, 

:  In  this  case  the  original  LqDV^  and  MqDV^  are  not  changed 

but  RqDV^^  is  modified  and  must  be  recalculated 0  The 
original  cross-section  appears  in  Figure  25 •  Using  formulas 
of  subroutines  RJ13  and  RJl4  the  new  volume  P  is  P  ~  volume 
FGED  +  volume  ADEC  -  volume  ABGF * 

:  In  this  case  the  original  RqDV^  is  zero*  the  original  LqDV^ 

is  unchanged  and  the  modified  MqDV^  is  computed  from  the 
formula  for  subroutine  RJ13« 

:  In  this  case  the  original  MqDV^  and  RqDVdu  are  zero  and  the 

modified  LqDV^  Is  computed  from  the  formulas  for  subroutines 

RJ13  and  14  as  the  difference  of  two  volumes* 

In  this  case  the  contribution  of  level  L,  on  the  u  th  search 

d 

turn  is  zero* 
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Figure  25:  Computation  of  modified  volume  RqDV^ »  Length  FG  is 

Bq0(d-1)  "  (u  “  1 ^  TES*  Area  BGH  represents  the  original 
volume  RqDVdu„  Line  AF  is  part  of  the  vertical  through 
the  weapon. 


125 


Note 

Note 

Note 

Note 

Note 

Note 

Note 

Note 


9:  The  contributions  of  levels  above  the  critical  depth  have 

now  been  determined* 

10s  The  contribution  of  level  is  not  changed* 

11 s  The  missed  volume  on  the  previous  search  turn  has  been  reduced 
leaving  less  to  be  picked  up  on  the  unchanged  second  turn* 

The  formulas  are  those  of  subroutines  RJ13  and  RJl4* 

12s  The  missed  volume  in  both  the  current  and  the  previous 

search  has  been  reduced*  The  formula  is  identical  to  that 
discussed  in  Note  7« 

13s  The  scanned  volumes  indexed  by  one  set  of  q#  u^  and  i  are 

now  computed*  The  next  step  is  to  change  u  and  compute  the 
next  set  of  volumes* 

lht  Scanned  volumes  on  the  i  th  ping  of  all  search  turns  above 
search  depth  are  now  computed*  The  next  step  is  to  compute 
the  volumes  below  search  depth* 

15:  The  increments  of  scanned  volume  for  ping  i ^  are  now 

computed*  The  next  step  is  to  compute  SV^* 

16s  Having  found  SViu>  the  index  is  changed  to  find 

SV,,  ,  x  . 

(i  +  l)u 
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Note  1:  The  problem  is  to  compute  areas  as  depicted  horizontally  and 
vertically  lined  in  Figure  l60  Using  Figure  16S  with  origin 
at  point  C  and  x  axis  toward  point  C,  the  equation  of  the 
target  volume  cylinder  circumference  on  ping  i  is 
(x  +  R  )2  +  y2  =  R.2 

O  1 

Changing  to  polar  coordinates  (r^  p)  with  the  same  origin 

and  measuring  angles  positive  counterclockwise  from  the 

x  axis,  the  equation  becomes 

r2  +  (2R  cosB  )  •  r=(R2  =R2)  -  o 

Solving  for  r, 

r  =  R.  ^  1  -  (R  /r.  )2  sin  2p  =R  cosP 
i  'o'  i  7  o 

To  find  the  area  for  ping  i*  r(r  d  p  )  is  integrated  from 
p'  +  (i  -  2)6  to  p*  +  (i  -  1)6  . 


Thus 

area  = 

Let  p*  +  (1-1)6  =  B 
p»  +  (1-2)6  =  A 


p'  + (1-1)6  . 

T  dp 

p*+ (1-2)6 


Then 
•B 
A 


f  r2  dp  =  f  [(R*  -  R02)  ♦  2  R  2  C0S2p 

-'a  A  L  °  ° - 5— -r- 

-  2  Rq^V  1  -  (Ro^i.)  8ln  p  C0SP 

JrB 

(Ri2  -  R02)  dp  =  (R±2  -  Ro2)(B  -  A) 

A 


dp 
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B 


2  R  2  cos2p  dp  =  R  2 
o  r  r  o 


p  +  (1/2)  sin2p 


A 


Let  x  -  sinp 
Then  dx  =  cosp  dp 


-2RQRiJ"  ^1  -  (^Ai)2  sin2p  cosp  dp 

=  -  2  R^  fB  Vl  -  (\A±)2  7  dx 

«<  A 


Let  sinp  =  x 

Then  cosp  dp  =  (Rq/R^)  dx 


-  2 


RoR 


jy 


(E0/S: 


n2  2 

)  x 


dx  =  -  2  R* 


r 

Jk 


2 

cos  p  d  p 


t 


p  +  (l/2)sin2p 


B 

A 


Therefore  * 


r2  dp  =  (Rq2  -  R±2)  (sin2B  -  sin2A)  /  2 
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ACQUISITION  PROBABILITIES:  I  PINGS 


x  ,  INPUT _ _  ...  INEPT — _  _ _ _ 

1  ) - (  APPENDIX  C7  OUTPUT  ) — (  Aw,  W«  ,  W»  ) - (1— »u 


\/ 

SVlu  PQi/TT^g^ 

- (  1  +  > — (i  : 

L(  1  U-g>U 


^ ~5>U  g(M-l)u^l//^)^  *  gl(u*l)fc  •“  j 


.u  ;  U}-^ - {  1  +  u  u  y~ 


1  ■*•  s  — (1  — >u 


O-^g 


MU 


HSlu^  ggj - (  l  *  S  — ( 1  + 1-»1 


1  +  P  ->U 


(W'/Aw)  -  (W«/am)  -5>  R 

- <1^r  >-| 

 1 

W'  — =>■  w. 


& 


[(wlr  “  wOSD^  PRVwlt 


I 


1->AP. 


[<f»s 


0H3>PA(''lr) 


>  £ - 

_  --  x  / 

Nip«^ip1 

— ( 1 4  ) — fi  :  ip 

1  -  APa  -5>  PA(«lr) 


— i 

OUTPUT 

1 

(  W’  wlr5  r  =  x-  R  ) 

■(I) 
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Note  Is  The  problem  constants  are  now  determined.  Next  the 


probability  P  is  determined  for 
Hjk 


j  =  1  >  °  °  °  >  ^i,  >  ^  =  1 , 


K. 
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